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Chapter 1 
Introduction 



It is well-known in classical probability theory that a expectation E satisfies 
the following relation for random variables X and Y: 

E[aX + Y + c] = aE[X] + E[Y] + c, Vo, c e M. 

A sublinear expectation E satisfies the following weaker condition: 

E[X + Y]< E[X] + E[Y], E[\a\X + c] = \a\E[X] + c. 

This E keeps the monotonicity property: If X > Y then E[X] > E[Y]. 

The notion of sublinear expectations is proved to be a basic tool in volatil- 
ity uncertainty which is crucial in superhedging, superpricing (See [3] and 
[44]) and measures of risk in finance which caused a great attention in fi- 
nance since the pioneer work of Artzner, Delbaen, Eber and Heath [1] and 
[2]. This is also the start point of a new theory stochastic calculus which 
gives us a new insight to characterize and calculate varies kinds of financial 
risk. 

In this note we will introduce a crucial notion of G-normal distributions 
corresponding to the well-known normal distributions in classical probabil- 
ity theory. This G-normal distribution will bring us to a new framework 
of stochastic calculus of Ito's type through the corresponding G-Brownian 
motion. We will also present analytical calculations and some new statis- 
tical methods with application to risk analysis in finance under volatility 
uncertainty. 

Our basic point of view is: sublinear expectation theory is very like its 
special situation of linear expectation in the classical probability theory. Un- 
der a sublinear expectation (or even more general nonlinear expectation) 
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space we still can introduce the notion of distributions, of random variables, 
as well as the notions of joint distributions, marginal distributions, etc. We 
still denoted X ~ Y if X and Y are identically distributed. We still have 
the notion of independence. A particularly interesting phenomenon in sub- 
linear situations is that a random variable Y is independent to X does not 
automatically implies that X is independent to Y. 

We will prove two important theorems in the framework of sublinear 
expectation theory: The law of large number and the central limit theo- 
rem. A very interesting result of our new central limit theorem under a 
sublinear expectation is that a sequence of zero-mean independent and iden- 
tically distributed random will converge in law to a 'G-normal distribution' 
A/"(0, [<l 2 ,<7 2 ]). Briefly speaking a random variable A in a sublinear expec- 
tation space is said to be A/"(0, [a 2 , cf 2 ]) -distributed if for any Y independent 
and identically distributed w.r.t. A and for any real function a and b we 
have aX + bY ~ Va 2 + b 2 X. Here a 2 = E[A 2 ] and a 2 = -E[-A 2 ]. In a 
special case a 2 = a 2 , this G-normal distribution becomes a classical normal 
distribution Af(0,a 2 ). 

We will define a sublinear expectation on the space of continuous paths 
from R + to R d which will be an analogue of Wiener's law, by which a 
G-Brownian motion is formulated. Briefly speaking a G-Brownian motion 
(B t )t>o is a continuous process with independent and stationary increments 
under a given sublinear expectation E[-]. 

G-Brownian motion has a very rich and interesting new structure which 
non trivially generalizes the classical one. We can establish the related 
stochastic calculus, especially G-Ito's integrals (see [35, 1942]) and the re- 
lated quadratic variation process (B). A very interesting new phenomenon 
of our G-Brownian motion is that its quadratic process (B) also has inde- 
pendent and stationary increments. The corresponding G-Ito's formula is 
obtained. We then introduce the notion of G-martingales and the related 
Jensen inequality for a new type of "G-convex" functions. We have also 
established the existence and uniqueness of solution to stochastic differential 
equation under our stochastic calculus by the same Picard iterations as in the 
classical situation. Books on stochastic calculus e.g., [Hj, [32], [31], [36], [4"T] . 
[IT] . [62] . [63] . [69] are recommended for understanding the present results 
and some further possible developments of this new stochastic calculus. 

A sublinear expectation can be regarded as a coherent risk measure. This, 
together with the related conditional expectations E[-|W t ]i>o makes a dy- 



namic risk measure: G-risk measure. 

The other motivation of our G-expectation is the notion of (nonlinear) 
g-expectations introduced in [50], [51]. Here g is the generating function of 
a backward stochastic differential equation (BSDE) on a given probability 
space (O, J-, P). The natural definition of the conditional (^-expectations 
with respect to the past induces rich properties of nonlinear g-martingale 
theory (see among others, [3], [7J, [8], [9], rj5], [16], rjO], [II], [3B], [39], [52], 
[56] . [57] . [59]). Recently (/-expectations are also studied as dynamic risk 
measures: g-risk measure (cf. [63], [5], [2]])- Fully nonlinear super- hedging 
is also a possible application (cf . [3] , [33] and [12] where new BSDE approach 
was introduced). 

The notion of g-expectation is defined on a given probability space. In 
[55] (see also [53])- As compared with the framework of ^-expectations, 
the theory of G-expectation is more intrinsic, a meaning similar to "intrinsic 
geometry" in the sense that it cannot be based on a given (linear) probability 
space. Since the classical Brownian expectation as well as many other linear 
and nonlinear expectations are dominated by our G-Expectation, our theory 
also provides a flexible theoretical framework. 

The whole results of this paper are based on the very basic knowledge 
of Banach space and the parabolic partial differential equation ( 13.31) . When 
this G-heat equation f l 3 . 3 H is linear, our G-Brownian motion becomes the 
classical Brownian motion. This paper still provides an analytic shortcut to 
reach the sophistic Ito's calculus. 

A very basic knowledge of Banach space and probability theory are neces- 
sary. Stochastic analysis, i.e., the theory of stochastic processes will be very 
helpful but not necessary. We also use some basic knowledge on the smooth 
solutions of parabolic partial differential equations. Good background in 
statistics will be very helpful. This note was written for several series of 
lectures: In the 2nd Workshop "Stochastic Equations and Related Topic" 
Jena, July 23-29, 2006; In graduate Courses of Yantai Summer School in Fi- 
nance, Yantai University, July 06-21, 2007; as well as in Graduate Courses of 
Wuhan Summer School, July 24-26, 2007. Also in mini-course of Institute of 
Applied Mathematics, AMSS, April 16-18 2007 and a mini-course in Fudan 
University, May 2007; In graduate courses of CSFI, Osaka University, May 
15- June 13, 2007. The hospitalities and encouragements of the above insti- 
tutions and the enthusiasm of the audiences are the main engine to realize 
this lecture notes. 

I thank for many comments and suggestions given during those courses, 
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especially to Li Juan and Hu Mingshang. References are given at the end. 
Historical remarks are still under preparation. This lecture note are mainly 
based on my recent research papers: 

Peng, S. (2006) G-Expectation, G-Brownian Motion and Related Stochastic 
Calculus of Ito's type, (pdf-file available in arXiv:math.PR/0601035v2, to 
appear in Proceedings of the 2005 Abel Symposium, Springer. 

Peng, S. (2006) Multi-dimensional G-Brownian motion and related stochastic 
calculus under G-expectation, Preprint, (pdf-file available in arXiv:math. PR/060 1699 
v2). 

Peng, S. (2007) Law of large numbers and central limit theorem under non- 



linear expectations, in arXiv:math.PR/0702358vl 13 Feb 2007 



which were stimulated by: 

Peng, S. (2004) Filtration Consistent Nonlinear Expectations and Evalua- 
tions of Contingent Claims, Acta Mathematicae Applicatae Sinica, English 
Series 20(2), 1-24. 

Peng, S. (2005) Nonlinear expectations and nonlinear Markov chains, Chin. 
Ann. Math. 26B(2) ,159-184. 



Chapter 2 

Risk Measures and Sublinear 
Expectations 

2.1 How to measure risk of financial positions 

Let Q be a set the set of scenarios. We are given a linear subspace X of real 
valued and bounded functions on Q X is the collection of are all possible risk 
positions in a financial market. We assume that all constants are in X and 
that X G X implies \X\ EX. For each X G X we denote 

X* = supX(V) and ||X|| = |X|*. 

||-|| is a Banach norm on X. In this section X is assumed to be a Banach 
space, i.e., if a sequence {X,{\°2 =1 of X converges uniformly to some function 
X on Q, then X E X. 

Remark 2.1.1 If S E X , then for each constant c, S\/ c, S Ac are all in X. 
One typical example in finance is that S is the tomorrow 's price of a stock. 
In this case any European call or put options of forms 

(S-k) + , {k-S) + 

are in X . 

Problem 2.1.2 Prove that if S G X then ip(S) G X for each continuous 
function (p onM. 

9 
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Example 2.1.3 Let (^,^ r ) be a measurable space and let L°°(J7, T) be the 
space of all bounded J- '-measurable random variables. L°°(r2,jF) is considered 
as a space of risk positions in a financial market. L°°(J7,jF) is a Banach 
space under the norm H-X^l^ = sup wgC |X(u;)|. 

Example 2.1.4 If Q is a metric space, then we can consider X = Cb(Q), 
the set of all bounded and continuous functions on Q. 

2.1.1 Coherent measures of risk 

A risk supervisor is responsible for taking a rule to tell traders, stocks com- 
panies, banks or other institutions under his supervision, which kind of risk 
positions is unacceptable and thus a minimum amount of risk capitals should 
be deposited to make the positions to be acceptable. The collection of ac- 
ceptable positions is defined by: 

A = {X G X, X is acceptable}. 

This set has the following economically meaningful properties: 

Definition 2.1.5 (Coherent acceptable set) 
(i) Monotonicity: 

X eA,Y>X => Y eA 

(ii) e A but -I £ A. 

(Hi) Positively homogeneity 

X eA => XX g A, VA > 0. 

(iv) Convexity: 

X, Y G A => aX + (1 - a)Y G A, a G [0, 1] 

Remark 2.1.6 (Hi) and (iv) imply 
(v) Sublinearity: 

X, Y e A =► fiX + vY e A, V>, v > 0. 
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Remark 2.1.7 If we remove the condition of the positive homogeneity, then 
A is called a convex acceptable set. In this course we mainly study the coher- 
ent case. Once the rule of the acceptable set is fixed, the minimum require- 
ment of risk deposit is then automatically determined. 

Definition 2.1.8 (Risk measure related to a given acceptable set 

A) The functional p(-) defined by 

p(X) = p A (X) := inf{m Gl: m + X e A}, X eX 
is called the risk measure related to A. 
It is easy to see that 

p(X + p(X)) = 0. 

Proposition 2.1.9 p(-) is a coherent risk measure, namely 

(a) Monotonicity: X > Y implies p(X) < p(Y); 

(b) Constant preservation: p(l) = — p{— 1) = — 1; 

(c) Sublinearity: For each X,Y G X, p(X + Y) < p{X) + p(Y); 

(d) Positive homogeneity: p(XX) = Xp(X), VA > 0. 

Proof, (a), (b) are obvious. To prove (c) we first have p(pX) = pp(X). In 
fact the case p = is trivial; when p > 0, 

Tfl 

p(fjtX) = inf {m 6l: m + pX e A} (n = — ) 

= a inf {n 6l: n + X G A} = pp(X). 

Now, for each positive a and is, 

p(pX + uY) = inf {m G R : m + (pX + uY) G A} 

= inf {m + n : m,n Gl, (m + pX) + (n + vY) G A} 

< inf{m eR: m + pX G A} + inf {n Gl: n + uY e A} 

= P {pX)+p{vY). 
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■ 

Now let a p(-) be a functional satisfying (a)-(d). Then we can inversely 
define 

A p = {X E X : p(X) < 0}. 

It is easy to prove that A p satisfies (i)-(iv). 

2.1.2 Sublinear expectation of risk loss 

From now we will denote X to be a loss position, namely —X is the corre- 
sponding financial position. Related to a coherent risk measure p, we evaluate 
the risk loss X by: 

E[X] := p(-X), X EX. 

This functional satisfies the following properties: 

(a) Monotonicity: 

X > Y => E[X] > E[Y]. 

(b) Constant preserving 

E[ c ] = c, Vc e R. 

(c) Sub-additivity: For each X, Y e X, 

E[X + Y] <E[X]+E[Y]. 

(d) Positive homogeneity: 

E[XX] = \E[X], VA > 0. 

A real valued functional E[-] defined on X satisfying (a)-(d) will be called 
a sublinear expectation and will be systematically studied in this lecture. 

Remark 2.1.10 (c) + (d) is called sublinearity. This sublinearity implies: 

(e) Convexity: 

E[aX + (1 - a)Y] < aE[X] + (1 - a)E[Y], Vet G [0, 1]; 
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Remark 2.1.11 (b) + (d) implies 
(f) Cash translatability: 

ELY + c] = E[X] + c. 

Indeed, we have 

E[X] +c = ELY] - E[-c] 

< ELY + c}< ELY] + E[c] = ELY] + c. 

Remark 2.1.12 (c) + (d) <^=^> (d) + the following Convexity: 

E[aX + (1 - a)Y] < aE[X] + (1 - a)E[Y], Va G [0, 1]. 

Remark 2.1.13 We recall the notion of the above expectations satisfying 
(c)-(d) was systematically introduced by Artzner, Delbaen, Eber and Heath 
|iy, JBi, in the case where Q is a finite set, and by Delbaen J2U^ in general 
situation with the notation of risk measure: p{X) = E[—X]. See also in 
Huber IS^tf for even early study of this notion E (called upper expectation E* 
in Ch.10 of [33]) in a finite set Q. See Rosazza Gianin [6$ or Peng f53 ], El 
Karoui & Barrieu f^/, [2^ for dynamic risk measures using g -expectations. 
Super-hedging and super pricing (see I25lj and JEEP are also closely related to 
this formulation. 

2.1.3 Examples of sublinear expectations 

Let E 1 and E 2 be two nonlinear expectations denned on (Q,X). E 1 is said 
to be dominated by E 2 if 

E X [X] -E 1 ^] <E 2 [X-F], \/X,YeX. 

The strongest sublinear expectation on X is 

E°°[X] :=X* = sup X(lu). 

Namely, all other sublinear expectations are dominated by E°°[-]. From (c) 
a sublinear expectation is dominated by itself. 

Notation 2.1.14 Let Vf be the collection of all finitely additive probability 
measures on (tt,J-). 
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Example 2.1.15 (A linear expectation) We consider U^(VL, J 7 ) the col- 
lection of risk positions with finite values. It is a subspace of X consisting of 
risk positions X of the form 

N 

X(u)= J^SiWw), x t eR, A i eJ 7 ,i = l,--- ,N. (2.1) 

i=i 

It is easy to check that, under the norm ||-|| , L^°(f2, J 7 ) is dense on L°°(f2, J 7 ). 
For a fixed Q G Vf and X G Lg°(0, T) we define 

N N „ 

E Q [X] = E Q ^T Xi l Ai {uj)\ := Y>Q(A) = / X{u)Q{dw) 
i=i i=i -^ 

£"g : L^°(fi,jF) -4 R is a linear functional. It is easy to check that Eg 
satisfies (a)-(b). It is also continuous under U^X" || 00 . 

\E Q [X]\< sup |X(w)| = iixil. 

Since Lg° is dense in L°° we £/ien can extend Eq from L§° to a linear con- 
tinuous functional on L°°(fi, J 7 ). 

Proposition 2.1.16 The functional Eq[-] :X \— ► R satisfies (a) and (b). 
Inversely each linear functional rj(-) :X \— ► R satisfying (a) and (b) induces 
a finitely additive probability measure via Q n {A) = 77(1,4), A G T . The 
corresponding expectation is 77 itself: 



r](X) = / X(w)Q„(du;). 

Notation 2.1.17 Lei Q EVf be given and let X : Q 1— ► R &e a J 7 -measurable 
function such that \X(u)\ < 00 /or eac/i c<j. T/je distribution of X on 
(fi, J 7 , Q) is defined as a linear functional 

F x [<p] = E Q [<p(X)] : if G L°°(R,8(R)) ^ R. 

Fx[-] is a linear functional satisfying (a) and (b). Thus it induces a finitely 
additive probability measure on (R, £>(R)) via Fx(B) := Fx(1b) , B G £>(R). 
We have 

F x [<p]= f <p(x)F x (dx). 

Jr 
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Remark 2.1.18 Usually, people call the probability measure Fx(-) to be the 
distribution of X under Q. But we will see that in sublinear or more general 
situation the functional version is necessary. 

Definition 2.1.19 Let E be a linear finitely additive probability defined on 
a measurable space (£l,J-). We say that a random variable Y is independent 
to X if, for each cp G L°°(M 2 ) 



E[<p(X,Y)]=E\E[<p(x,Y)] 



x=X\- 



2.1 A Representation of a sublinear expectation 

Let a risk supervisor take a linear expectation E[-] to be his risk measure. 
This means that he takes the corresponding finitely additive probability in- 
duced by E[-} as his probability measure. But in many cases, he cannot 
decide precisely which probability he should take. He has a set of finitely 
additive probabilities Q. The size of Q characterizes his model-uncertainty 
A robust risk measure under such model uncertainty is: 

E S [X] = sup E Q [X] : X h-> R. 
QeQ 

It is easy to prove that E s [-] is a sublinear expectation. 

Theorem 2.1.20 A sublinear expectation E[-] has the following representa- 
tion: there exists a subset Q cP/, such that 

ELY] = max E Q [X], VX G X. 

Proof. It suffices to prove that, for each Xo G X, there exists a Qx G Vf 
such that E Qxq [X] < E[X], for all X, and such that Q Xo (X ) = E[X ]. 

We only consider the case E[X ] = 1 (otherwise we may consider X = 
X -E[X ] + l). 

Let U\ := {X : E[X] < 1}. Since U\ is an open and convex set, by 
the well-known separation theorem for convex sets, there exists a continuous 
linear functional rj on X such that i](X) < i](Xq), for all X G U\. Since 
G U\ we have particularly i](Xq) > 0, and we then can normalize rj so that 
r)(X ) = 1 = ELYo]. Thus rj satisfies: 

r](X) < 1, VX G X, s.t. E[X] < 1. 
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By the positive homogeneity of E, for each c > 0, 

7](X) < c, VI6 X, s.t. E[X] < c. 

This implies that 

E[X] > 7]{X), for all X such that E[X] > 0. (2.2) 

We now prove that r] is monotone in the sense of (a). For a given Y > 
and A > 0, since — XY G U\, we have 

V(-XY) < 7](X ), thus T](Y) > -X^riiXo) = -A -1 , VA > 0. 

From which it follows that 77(F) > and thus (a) holds true. From (12.21) we 
have 77(1) < E[l] = 1. We now prove that 77(1) = 1. Indeed, for each c > 1, 
we have 

E[2X - c] = 2E[X ] -c = 2-c<l, 

hence 2X — c & U\ and 

77(2X - c) = 2 - C7](l) < 1 orn(l)>-, Vol; 

c 

hence 77(1) = 1. This, together with (12.21) . yields 

E[X] > 7](X), for all X EX. 
We thus proved the desired result. 



Chapter 3 

LLN and Central Limit 
Theorem 



3.1 Preliminary 



The law of normalsize numbers (LLN) and central limit theorem (CLT) are 
long and widely been known as two fundamental results in the theory of 
probability and statistics. A striking consequence of CLT is that accumu- 
lated independent and identically distributed random variables tends to a 
normal distributed random variable whatever is the original distribution. It 
is a very useful tool in finance since many typical financial positions are ac- 
cumulations of a large number of small and independent risk positions. But 
CLT only holds in cases of model certainty. In this section we are inter- 
ested in CLT with variance-uncertainty We will prove that the accumulated 
risk positions can converge 'in law' to what we call G-normal distribution, 
which is a distribution under sublinear expectation. In a special case where 
the variance-uncertainty becomes zero, the G-normal distribution becomes 
the classical normal distribution. Technically we introduce a new method to 
prove a CLT under a sublinear expectation space. 

In the following two chapters we will consider the following type of spaces 
of sublinear expectations: Let Q be a given set and let Ti be a linear 
space of real functions defined on Q such that if X\, ■ ■ ■ ,X n G Ti then 
tp(Xi,--- ,X n ) G Ti for each (p G Ci.Li P (^ n ) where Ci.Li P {^ n ) denotes the 

17 
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linear space of functions tp satisfying 

\<p(x) - v(y)\ < C(l + \x\ m + \y\ m )\x - y\, \fx,y e R n , 
for some C > 0, m G N depending on ip. 

Ti is considered as a space of "random variables" . 

Remark 3.1.1 In particular, if X,Y G Ti, then \X\, X m G Ti are in Ti. 

More generally <p(X)if)(Y) is still in Ti if (p,ip G Cj.£* p (R). 

Here we use Ciupi^ 1 ) m our framework only for some convenience of 
techniques. In fact our essential requirement is that Ti contains all constants 
and, moreover, X G 7i implies \X\ EH. In general Ci.Li P (M n ) can be replaced 
by the following spaces of functions defined on ¥L n . 

• L°°(]R ra ): the space bounded Borel-measurable functions; 

• Cb{M. n ): the space of bounded and continuous functions; 

• C5 (W a ): the space of bounded and fc-time continuously differentiate 
functions with bounded derivatives of all orders less than or equal to 
k; 

• C un if(M. n ): the space of bounded and uniformly continuous functions; 

• Cb.LipiJ^ 1 )'- the space of bounded and Lipschitz continuous functions; 

• L (R n ): the space of Borel measurable functions. 

Definition 3.1.2 Sublinear expectation^ on Ti is a functional E : Ti i— > 
R satisfying the following properties: for all X, Y G Ti, we have 

(a) Monotonicity: IfX>Y then E[X] >E[Y}. 

(b) Constant preserving: E[c] = c. 

(c) Sub-additivity: ' E[X] - E[Y] < E[X - Y\. 

(d) Positive homogeneity: E[XX\ = \E[X], VA > 0. 

(In many situation (c) is also called property of self-domination) . The triple 
(fl, Ti, E) is called a sublinear expectation space (compare with a proba- 
bility space (Q, T , P)^. 
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Example 3.1.3 In a game we select at random a ball from a box containing 
W white, B black and Y yellow balls. The owner of the box, who is the banker 
of the game, does not tell us the exact numbers ofW,B and Y. He or she 
only informs us W + B + Y = 100 and W = B G [20, 25]. Let £ be a random 

variable 

{1 if we get a white ball; 

if we get a yellow ball; 

—1 if we get a black ball. 

Problem: How to measure a loss X = </?(£), for a given function tp on R. We 
know that the distribution of £ is 



p i m p 



1 1 I . , . r 2 -2n 

2 1 ~P 2 



with uncertainty: p G[cx 2 ,(T ] = [0.4,0.5]. 



Thus the robust expectation of X = </?(£) is: 
%(0]:=sup£ P b(0] 



Pev 



= sup [^(l)+^(-l)] + (l-pV(0)]. 

£ /ias distribution uncertainty. 

Example 3.1.4 A more general situation is that the banker of a game can 
choose among a set of distribution {F(0, ^4)}AeB(M),0ee of a random variable 
£ . In this situation the robust expectation of a risk position </?(£) for some 

¥ e C LLip (R) is: 

1^(0] := sup f <p(x)F(6,dx). 

3.2 Distributions and independence 

We now consider the notion of the distributions of random variables under 
sublinear expectations. Let X = (Xi, ■ ■ ■ ,X n ) be a given n-dimensional 
random vector on a sublinear expectation space (fii,7ii,E). We define a 
functional on C/.Li P (R n ) by 

¥ x [ip] := E[<p(X)] : y? G C LLip (R n ) h-> (-oo, oo). (3.1) 

The triple (R n , C;.Lj P (R n ),Fx[-]) forms a sublinear expectation space. Fx is 
called the distribution of X. 
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Definition 3.2.1 Let X\ and X2 be two n-dimensional random vectors de- 
fined respectively in sublinear expectation spaces (Qi, Tii, Ei) and (Q2, ^2, E2). 
They are called identically distributed, denoted by X\ ~ X2, if 

E 1 [ l p(X 1 )} = E 2 [<p(X 2 )}, Vy? G C LLip (R n ). 

It is clear that X\ ~ X 2 if and only if their distributions coincide. 

Remark 3.2.2 If the distribution Fx of X G 7i is not a linear expectation, 
then X is said to have distributional uncertainty. The distribution of X has 
the following four typical parameters: 

■p:=E[X], fi:=-E[-X], a 2 := t[X 2 ], a 2 := -t[-X 2 ]. 

The subsets [/i,/l] and [a 2 ,? 2 ] characterize the mean-uncertainty and the 
variance-uncertainty of X . The problem of mean uncertainty have been stud- 
ied in [Chen-Epstein] using the notion of g- expectations. In this lecture we 
are mainly concentrated on the situation of variance-uncertainty and thus set 

JI = j2. 

The following simple property is very useful in our sublinear analysis. 

Proposition 3.2.3 Let X,Y e H be such that E[Y] = -E[-Y], i.e. Y has 
no mean uncertainty. Then we have 

E[X + Y] =E[X]+E[Y]. 

In particular, ifE[Y] = E[-Y] = 0, then E[X + Y}= E[X]. 

Proof. It is simply because we have E[X + Y]< E[X] + E[Y] and 

E[X + Y]> E[X] - E[-Y] = E[X] + E[Y]. 

■ 
The following notion of independence plays a key role: 

Definition 3.2.4 In a sublinear expectation space (Q, 7i, E) a random vector 
Y = (Y~i, ■ ■ ■ ,Y n ), Yi G Ti is said to be independent to another random 
vector X = (Xi, ■ ■ ■ ,X m ), Xi G 7i under E[-\ if for each test function ip G 
C LLi p(R m x R n ) we have 

&[<p(X,Y)]=±\E[<p(x,Y)] x=x ]. 
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A random variable Y G 7i is said to be weakly independent to X G Ji® m if 
the above test functions tp are only taken from the following class: 

v {x,y)=Mx)+Mx)y + Mx)y 2 , A e c LLip (w n x r), 2 = 0,1,2. 

Remark 3.2.5 In the case of linear expectation, this notion of independence 
is just the classical one. It is important to note that under sublinear expecta- 
tions the condition l Y is independent to X " does not implies automatically 
that "X is independent to Y". 

Example 3.2.6 We consider a case where X, Y G ri are identically dis- 
tributed and E[X] = E[-X] = but a 2 = E[X 2 ] > a 2 = -E[-X 2 ]. We 
also assume that E[\X\] = E[X + + X~] > 0, thus E[X + ] = |E[|X| + 
X] =|E[|X|] > 0. In the case where Y is independent to X , we have 

E[XY 2 ] = E[X + a 2 - X-a 2 ] = (a 2 - a 2 )E[X + ] > 0. 
But if X is independent to Y we have 

E[XY 2 ] = 0. 

The independence property of two random vectors X, Y involves only the 
joint distribution of (X,Y). The following result tell us how to construct 
random vectors with given sublinear distributions and with joint indepen- 
dence. 

Proposition 3.2.7 Let Xi be Hi- dimensional random vectors respectively in 
sublinear expectation spaces (Qi,rii, Ej), % = 1, • • • ,N. We can construct 
random vectors Y"i, • • • , Y^ in a new sublinear expectation space (Q, 7Y,E) 
such that Yi ~ Xi and such that Y i+i is independent to (Y 1; ■ • • ,Yj), for each 
i. 

Proof. We first consider the case N = 2. We set: tt = R" 1 x IR™ 2 , i.e., 
uj = {x\ ) X2) G Q, X\ G IR™ 1 X2 G M™ 2 ; a space of random variables 7i = 
{X(u) = tp(uj),tp G CixipiS^)} and a functional E on (fi,TC) defined by 



E[^(y)]=Ei[^(Xi)], where <p( Xl ) := E 2 [i P (x 1 ,X 2 )}, Xl G 

VF(cu) = <p(w), ^ G C,L ip (K ni x R" 2 ). 



pni 



• 
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It is easy to check that E forms a sublinear expectation on (fi,7i). Now let 
us consider two random vectors in 7i: 



Yl(u;) = xi, F 2 (cj) = x 2 , w = (xi, x 2 ) 6 M ni x 



?"2 



It is easy to check that Y±, Y 2 meet our requirement: Y 2 is independent to 
Y\ under E, and Y\ ~ Xi, Y 2 ~ X 2 . The case of A^ > 2 can be proved by 
repeating the above procedure. ■ 

Example 3.2.8 We consider a situation where two random variables X and 
Y in 7i are identically distributed and their common distribution is 



W x [<p] = ¥ Y [ip] = sup / <p(y)F(9,dy), <p 6 C LLip (R). 
see Jr 

where, for each 9 G 0, {F{6, A)}a&b(r) is a probability measure on (R, B(M,)). 
In this case l Y is independent to X " means that the joint distribution of X 
and Y is: 



^x,yW\ = SU P 



i> e C LUp (R 2 ) 



sup / ip(x,y)F(9 2 ,dy) 



F($ 1 ,dx), 



Remark 3.2.9 The situation l Y is independent to X" often appears when 
Y occurs after X , thus a very robust expectation should take the information 
of X into account. 

Definition 3.2.10 A sequence of random variables {^}^ 1 in H is said to 
converge in distribution under E if for each bounded and ip G Cb.up( 



oo 

converges. 



i=i 



3.3 G-normal distributions 

3.3.1 G-normal distribution 

A fundamentally important distribution in sublinear expectation theory is 
A/(0; [a 2 , a 2 ])-distributed random variable X under E[-]: 
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Definition 3.3.1 (G -normal distribution) In a sublinear expectation space 
(p,,E,TC), a random variable X G Ti with 

a 2 = E[X 2 ], a 2 = -E[-X 2 ] 

is said to be A/"(0; [a 2 , a 2 }) -distributed, denoted by X ~ A/"(0; [g^,a 2 ]), if for 
each Y G Ti which is independent to X such that Y ~ X we have 

aX + bY ~ Va 2 + b 2 X, Vo, 6 > 0. (3.2) 

Proposition 13.2.71 tells us how to construct the above Y based on the 
distribution of X. 

Remark 3.3.2 By the above definition, we have y/2E[X] = E[X + Y] = 
2E[X] and \/M[-X] = E[-X -Y] = 2E[-X] it follows that 

E[X] = E[-X] = 0, 

Namely an A/"(0; [a 2 , a 2 ]) -distributed random variable X has no mean uncer- 
tainty. 

Remark 3.3.3 If X is independent to Y and X ~ Y, such that A3.2\) sat- 
isfies, then —X is independent to —Y, —X ~ —Y. We also have a(—X) + 
b(-Y) ~ Va 2 + b 2 (-X), a,b>0. Thus 

X~M(0;\a 2 ,a 2 }) iff -X^Af(0;[a 2 ,W 2 }). 

The following proposition and corollary show that A/"(0; [a 2 , a 2 ]) is a 
uniquely defined sublinear distribution on (M, Cm p (M.)). We will show that 
an A/"(0; [a 2 , a 2 ]) distribution is characterized, or generated, by the following 
parabolic PDE defined on [0, oo) x R: 

d t u - G(dl x u) = 0, (3.3) 

with Cauchy condition u\ t =o = (p, where G, called the generating function 
of (I3.3p . is the following sublinear real function parameterized by a and a: 

G{a) := -E[X 2 a] = -(a 2 a + - a 2 «-), a G K. 

Here we denote a + := max{0, a} and a - := (— a) + . ( 13. 3D is called generating 
heat equation, or G-heat equation of the sublinear distribution A/"(0; [a 2 , a 2 ]). 
We also call this sublinear distribution G-normal distribution. 
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Remark 3.3.4 We will use the notion of viscosity solutions to the generating 
heat equation A3. 3]) . This notion was introduced by Crandall and Lions. For 
the existence and uniqueness of solutions and related very rich references we 
refer to Crandall, Ishii and Lions [11]. We note that, in the situation where 
a 2 > 0, the viscosity solution $3.3]) becomes a classical C 1+i % ,2+a -solution (see 
LJEjj , LJ3^ and the recent works of J^ and ffllj). Readers can understand 113.3]) 
in the classical meaning. 

Definition 3.3.5 A real-valued continuous function u G C([0, T] x R) is 
called a viscosity subsolution (respectively, supers olution) of $3.3]) if, for each 
function ip G C\ ((0, oo) x R) and for each minimum (respectively, maximum) 
point (t, x) G (0, oo) x R of ip — u, we have 

d t ip — G(dl x ip) < (respectively, > 0). 

u is called a viscosity solution of $3.3]) if it is both super and subsolution. 

Proposition 3.3.6 Let X be an A/"(0; [<l 2 ,cx 2 ]) distributed random variable. 
For each ip G Cj.z,ip(R) we define a function 

u(t,x) :=t[ip(x + VtX)], (t,x) G [0, oo) x R. 

Then we have 

u(t + s,x)='E[u(t,x + y/sX)], s>0. (3.4) 

We also have the estimates: For each T > there exist constants C, k > 
such that, for all t,s G [0, T] and x, y G R, 

\u{t,x)-u{t,y)\ < C(l + \x\ k + \y\ k )\x-y\ (3.5) 

and 

\u{t,x)-u(t + s,x)\ <C(l + \x\ h )\s\ 1/2 . (3.6) 

Moreover, u is the unique viscosity solution, continuous in the sense of A3. 5]) 
and IT3~o)) . of the generating PDE IT373)) . 

Proof. Since 

u{t, x) - u{t, y) = E[(p(x + y/tX)] - E[ip(y + y/tX)] 
< E[ip(x + VtX) - ip(y + y/iX)] 

<E[Ci(l + |X| fc + |a;| fc + |3/| fc )|a;-3/|] 

<C(l + \x\ k + \y\ k )\x-y\. 
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We then have ( 13.51) . Let Y be independent to X such that X ~ Y. Since X 
is 7V(0; [cr 2 ,a 2 ])-distributed, then 

u(t + s,x) = E[ip(x + y/t + sX)] 

= E[ip(x + yfsX + VtY)} 

= E[E[<p(x + V7sz + VtY)] z=x ] 

= R[u(t,X+y/sX)]. 

We thus obtain (13.41) . From this and (I3.5P it follows that 

u(t + s, x) — u(t, x) = E[u(t, x + \fsX) — u(t, x)] 

<E[Ci(l + |a;| fc + |X| fc )|s| 1/2 |X|]. 

Thus we obtain ([33]). Now, for a fixed (£, x) E (0, oo) xl, let ip E C 6 1,3 ([0, oo) x 
R) be such that ip > u and tp(t,x) = u(t,x). By (13. 4p it follows that, for 
5 E (0,t) 

< E[ip(t -5,x + V6X) - ip(t, x)] 
< -d t ip(t, x)S + E[d x ip(t, x)VSX + -dl x ip{t, x)5X 2 ] + C5 3/2 

= -d t ^(t, x)S + E[^d 2 xx ij(t, x)5X 2 ] + C6 3/2 
= -d t ^(t, x)S + 6G(d 2 xx ij(t } x)) + C5 3/2 . 

From which it is easy to check that 

[d t ^-G(d 2 x M(t,x)<0. 

It follows that u is a viscosity supersolution of (13.31) . Similarly we can prove 
that u is a viscosity subsolution of (13.31) . ■ 

Corollary 3.3.7 If both X and X are 7V(0; [a 2 , a 2 }) distributed. Then X ~ 
X . In particular, X ~ —X. 

Proof. For each tp E Qxip(R) we set 
u(t,x) := E[tp(x + VtX)], u(t,x) := E[cp(x + VtX)}, (t,x) E [0, oo) x R. 
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By the above Proposition, both u and u are viscosity solutions of the G-heat 
equation (13.31) with Cauchy condition u\ t =o = u\t=o = f- It follows from the 
uniqueness of the viscosity solution that u = u. In particular 

&[<p(X)]=fc[<p(X)]. 

Thus X ~ X. ■ 

Corollary 3.3.8 In the case where a 2 = a 2 > 0, A/"(0; [cr 2 ,a 2 ]) is just the 
classical normal distribution J\f(0;a 2 ). 

Proof. In fact the solution of the generating PDE (13.31) becomes a classical 
heat equation 



dtu = Y d lx u > u \t=o = V 



where the solution is 



u ( t ' x ) = ^7^^= </?u/) ex p( — 7p^—) d y 

V27CO-H J-oo Zo ~ t 

thus, for each (p, 

1 f°° v 2 

E[<p(X)] = u(l, 0) = -=== / <p(y) exp(-^)dy. 

V lTXO~ J-oo Zcr 

■ 

In two typical situations the calculation of Ek/>(X)] is very easy: 

• (i) For each convex ip, we have 

1 f°° v 2 

E[<p(X)] = -== / <p(y) exp(-^)dy 

V27TO- Z J-oo l0 ~ 

Indeed, for each fixed t > 0, it is easy to check that the function 

u(t, x) := Ef^o; + y/tX)] is convex: 

u(t, ax + (1 — a)y) = E[(p(ax + (1 — a)y + VtX)] 

< at[ip(x + VtX)} + (1 - a)E[if(x + VtX)} 
= au(t, x) + (1 — a)u(t, y) 

It follows that (d 2 x u)~ = and thus the G-heat equation ( 13.31) becomes 

a W \z i 

o t u = —O xx u, u\ t=0 = (p. 
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• (ii) But for each concave if, we have, 

E[p(X)] = - 7 L= r <p(y)eM-£)dy 

In particular 

E[X] = E[-X] = 0, E[X 2 ] = a 2 , -E[-X 2 ] = a 2 

and 

E[X 4 ] = 6a 4 , -E[-X A ] = 6a 4 

3.3.2 Construction of G-normal distributed random vari- 
ables 

To construct an 7V(0; [a 2 , a 2 ])-distributed random variable £, let u = u^ be 
the unique viscosity solution of the G-heat equation (13. 3p with u v \t=a = if. 
Then we take Q = R 2 , H = C LLip (R 2 ), u = fay) G R 2 . The corresponding 
sublinear expectation E[-] is defined by, for each X{uo) = ip{uj) such that 
i> e C LLip (R 2 ), 

E[X] = u^W(l,0), where we set $fa) := «^ (a,,0 (l,0). 

We now consider two random variables £(u;) = x, r/(u;) = y. It is clear that 
the 

E[<p(£)] = E[<p(r})] = v?{l, 0), Vipe C LLip {R). 

By the definition rj is independent to £ and £ ~ 77 under E. To prove that the 
distribution of £ and 77 satisfies condition ( 13.21) . it suffices to observe that, 
for each if e Ci,Li P (R) and for each fixed A > 0, x G R, since the function t> 
defined by t> (£, x) := u^^Xt, x + yXx) solves exactly the same G-heat equation 
(13.31) but with Cauchy condition v\ t=0 = if(x + vA x •), we then have 

E[<f(x + V\0] = v(l, x) = u v(VXx °(l, x) = u v (X, x), x£l. 
Thus, for each t > and s > 0, 

E[<p(y/i£ + y/sq)] = E[E[(f(Vix + V^) W 

= u u ^ s '-\t,0) = u v {t + s,0) 

= E[f{VtT^0]- 
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Namely >/££ + */sr] ~ \Jt + s£. Thus £ and r\ are both A/"(0; [cr 2 ,^? 2 ]) dis- 
tributed. 

We need to check that the functional E[-] : Ci.LipiM) \— > M forms a sub- 
linear expectation, i.e., (a)-(d) of Definition 13. 1.21 are satisfied. Indeed, (a) is 
simply the consequence of comparison theorem, or the maximum principle of 
viscosity solution (see Appendix). It is also easy to check that, when if = c, 
then the unique solution of (13.31) is also u = c; hence (b) holds true, (d) 
also holds since u Xv = Xu^, A > 0. The sub-additivity (c) will be proved in 
Appendix. 

3.4 Central Limit Theorem 

Our main result is: 

Theorem 3.4.1 (Central Limit Theorem) Let a sequence {Xi] c ^ =l in 7i be 
identically distributed with each others. We also assume that, each X n+ \ is 
independent (or weakly independent) to (Xx, • • • , X n ) for n = 1, 2, • • • . We 
assume furthermore that 

E[X X ] = E[-Xi] = ; E[X 2 ] = a 2 , -E[-X 2 ] = a 2 , 

for some fixed numbers < a < a < oo. Then the sequence {S n / y/u}°^ =1 , of 
the sum S n — X\ + \- X n , converges in law to A/"(0; [a 2 , a 2 ]): 

lim E[<p(^=)] = E[^(0], V^ G lip b (R), (3.7) 

n->oo x /n 



where £ ~ JV(0; [a 2 , a 2 }). 

Proof. For a small but fixed h > 0, let V be the unique viscosity solution of 
d t V + G(d 2 xx V) = 0, (£, x) e[0,l + h]x R, V\ t=1+h = <p. (3.8) 

We have, according to the definition of G-normal distribution 



V(t, x) = E[p{x + y/l + h-t£)]. 

Particularly, 

V(h,0) = E[<p(£)], V(l + h,x) = <p(x). (3.9) 
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Since ( 13.81) is a uniformly parabolic PDE and G is a convex function, thus, 
by the interior regularity of 7 (see Wang [67], Theorem 4.13) we have 

IMIci+«/ 2 .a+«([o,i]xR) < oo, for some a E (0,1). 

We set S = ± and S = 0. Then 

n u 

ra-1 



7(1, v^5 n ) - 7(0, 0) = $^{7((z + 1)5, V6S i+1 ) - V(i8, V6S,)} 

i=0 

n-1 

= J2{^ V ^ + ^ ^^+1) " V ^' v^+i)] + [7(*<J, VlS i+1 ) - 7(i5, v 7 ^ 

i=0 

n-1 



i=0 

with, by Taylor's expansion, 



Jl = d t V(i6, V5SJ6 + ^ xx V(iS, V6S t )Xf +1 6 + d x V(iS, VsSi)x i+1 Vs 

11= I [d t V((i + p)6,V6S i+1 )-d t V(i6,VSS i+1 )]dp5 
Jo 

+ [d t V(i5, VSS i+1 ) - d t V{i8, V6Si)]S 

+ f f [# B m V~5S t , + yfiX^y/S) - d 2 xx V(i8, y/6Si)YidpdrtXZv6. 
Jo Jo 



'0 Jo 
Thus 

n— 1 n— 1 n— 1 n— 1 

E[J2 4] - E[- £ 4] < E[7(l, vfa„)] - 7(0, 0) < E[j] 4] + £[£ 4] 

i=0 i=0 i=0 i=0 

(3-10) 
We now prove that E[^™ =0 ^s\ = 0- Indeed, the 3rd term of J\ has mean- 
certainty: 

E[d x V(iS, V6Si)X i+1 V6] = E[-d x V(iS, VdSi)X l+1 V5} = 0. 

For the second term, we have 

n\dl x V{i5, vfo 4 )X?u<$] = E[G(C7(z5, y/6Si))S\. 
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We then combine the above two equalities with d t V + G(d 2 x V) = as well 
as the independence of X i+ i to (Xi, • • • , X»), it follows that 

n— 1 n— 2 



eE^] = eE41 = - = o. 



i=0 i=0 

Thus (I3.10p can be rewritten as 

n— 1 n— 1 

',vl 



-E[- £ ij] < E[7(l, v^S n )] - 7(0, 0) < E[J] II 



i=0 i=0 

But since both d t 7 and 9^1^ are uniformly a-holder continuous in x and ~- 
holder continuous in t on [0, 1] x JR., we then have |/]| < C5 1+a ^ 2 [l + \X i+1 \ + 
\X i+1 \ 2+a }. It follows that 

E[|/j|] < C5 l+a/2 (l + E[|X x | a ] + E[|X!| 2+a ]). 

Thus 

_C(i)«/ 2 (l + E[|Xi| a + |Xi| 2+a ]) < E[7(l, v^5 n )] - 7(0, 0) 
n. 

< C(-) a/2 (l + ENX^ + |Xi| 2+a ]). 
n 

As n — > oo we have 

lim E[7(l, v^5 n )] = 7(0, 0). (3.11) 

On the other hand, for each t, t' G [0,1 + h] and x E R, 

|7(t,x)-7(t',x)| = |E[|^(x + Vl + h - tO - (p(Vl + h- t'0\] 



< k v \y/l + h-t-y/l + h-f\ x E[|£|] 

<cVl*-*'l, 

where fc v denotes the Lipschitz constant of <£>. Thus |7(0, 0) — 7(/i, 0)| < Cyh 
and, by (J3U), 

Ie^V^-e^V^s;)]! 

= |E[7(1, V6S n )} - E[7(l + fc, v^5 n )]| < Cv 7 ^. 
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It follows form lETTTj) and (GLHD that 

limsup |E[p(-p)] - Efo>(£)]| < 2Cv^. 



Since /z can be arbitrarily small we thus have 

limE[^)] =%>(£)]. 
in 



Corollary 3.4.2 The convergence \3.7\) holds for the case where (p is a 
bounded and uniformly continuous function. 

Proof. We can find a sequence {<fk}kLi in lipi,(M) such that (pk — > v 9 uni- 
formly on K. By 

l%(%] - Eb(0]l < ft?(%] - n<Pk(-%)]\ 






We can easily check that (13.71) holds. 
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Chapter 4 

G-Brownian Motion: 
1-Dimensional Case 

4.1 1- dimensional G-Brownian motion 

In this lecture I will introduce the notion of G-Brownian motion related to 
the G-normal distribution in a space of a sublinear expectation. We first give 
the definition of the G-Brownian motion. 

Definition 4.1.1 A process {Bt(ui)}t>o in a sublinear expectation space (Q, Ti, E) 
is called a G-Brownian motion if for each n G N and < ti, • • ■ ,t n < 

oo ; B tl , • • • , B tn G Ti and the following properties are satisfied: 

(i)B o (u) = 0; 

(ii) For each t, s > 0, the increment B t+S — B t is A/"(0; [a 2 s, a 2 s])- distributed 

and is independent to (B tl , B t2 , ■ ■ ■ , B tn ), for each n G N and < t± < • • ■ < 

t n <t. 

Remark 4.1.2 The letter G indicates that the process B is characterized by 
it's 'generating function' G defined by 

G(a) := E[aBf], aGl 

We can prove that, for each A > 0, (\~2B\ t )t>o is also a G-Brownian motion. 
For each to > 0. (B t+to — B to )t>o is also a G-Brownian motion. This is the 
scaling property of G-Brownian motion, which is the same as that of the usual 
Brownian motion. 

In this course we assume without loss of generality that a = 1 and q_ = 
a< 1. 

33 
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Theorem 4.1.3 Let (B t )t>o be a process defined in a sublinear expectation 
space (Q, Ti, E) such that 
(i)B o = 0; 

(ii) For each t, s > 0, the difference B t+S — B t and B s are identically dis- 
tributed and independent to (B tl , B t2 , • • • , B tn ), for each n GN and < t± < 

■■■ <t n <t. 

(iii) B = 0, E[B t ] =E[-B t ] = and lim ti0 E[| J B t | 3 ]t' 1 = 0. 

Then B is a Ga^-Brownian motion with a 2 —E[Bf] and a 2 = — E[— B 2 ]. 

Proof. We need only to prove that B t is A/"(0; [cr 2 £, ^^-distributed. We 
first prove that 

E[E 2 ] = a 2 t and - E[-B 2 ] = a 2 t. 

We set b(t) := E[B 2 }. Then 6(0) = and lim no b(t) < E[\B t \ 3 }^ -► 0. Since 
for each t, s > 0, 

b(t + s)= E[(B t+s ) 2 } = E[(B t+s -B s + B s ) 2 } = E[(B t+s - B s ) 2 + B 2 } 
= b(t)+b(s). 

Thus b(t) is linear and uniformly continuous in t; hence there exists a constant 
a > such that E[5 2 ] = a 2 t. Similarly, there exists a constant g_ G [0,<r] 
such that -E[-B 2 ] = a 2 t. We have a 2 = E[B 2 } > -E[-B 2 ] =a 2 . 

We now prove that B t is A/"(0; [a 2 £, <r 2 £])-distributed. We just need to 
prove that, for each fixed ip G Q.lj p (IR), the function 

u(t, x) := E[<p(x + B t )], (t, x) G [0, oo) x R 

is the viscosity solution of the following G>,— heat equation 

d t u - G(d 2 xx u) = 0, for t > 0; u\ t=0 = (p. (4.1) 

withG(a) = Gg-^(a) = \{a 2 a + — a 2 a~). Thus for each A/"(0; [a 2 , a 2 ])-distributed 
£, we have 

E[<p(B t )]=E[<p(y/it)], \f(peQ. Up (R). 

Thus B t ~ Vti ~ W(0; [a 2 t, a 2 t}). 

We first prove that u is locally Lipschitz in x and locally ^-Holder in t. 
In fact, for each fixed t, u(t, •) £Ci.Li P (R) since 

\E[ip(x + B t )] - E[<p(y + B t )\\ < ®[\<p(x + B t ) - <p(y + B t )\] 

<\E[C(l + \x\ m + \y\ m + \B t \ m )\x-y\] 
<Ci(l + |a;| m + |j/| m )|a;-j/|. 
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For each 5 G [0, t], since B t — B$ is independent to Bg, we also have 

u(t,x)=E[<p(x + B s + (B t -B s )] 

= E[E[^y + (B t -B 5 ))] y=x+Ss ], 



hence 
Thus 



u{t,x)=E[u{t-8,x + B s )}. (4.2) 

\u(t, x) — u(t — 8, x) | = |E[u(t — 5, x + B$) — u(t — 5, x)] \ 

< E[\u{t -8,x + B 5 )- u(t - 5, x) |] 
<E[C(l + \x\ m +\B 5 \ m )\B 5 \] 

< d{l + \x\ m )a5^ . 



To prove that u is a viscosity solution of ( 14.1ft . we fix a time-space (£, x) G 
(0, oo) x R and let v G C b ' ([0, oo) x R) be such that v > u and v(t,x) = 
u(t,x). From ( 14.21) we have 



v(t, x) = E[u{t -5,x + B s )\ < E[v{t -8,x + B s )] 

Therefore by Taylor's expansion, 

< E[v{t -8,x + B s )- v{t, x)} 
= E[v{t -8,x + B s )- v{t, x + B s ) + (v(t, x + B 5 ) - v{t, x))} 

= E[-d t v{t, x)8 + d x v{t, x)B 5 + ^d 2 xx v(t, x)B 2 5 + I s ] 

< -d t v(t, x)8 + E[-d 2 xx v{t, X )BJ\ - Ws] 
= -d t v(t,x)5 + G(d 2 xx (t,x))5-E[I 5 }. 
where 



-[d t v(t -{38,x + B s ) - d t v(t, x)]8d(5 
o 

+ / / [<%x v (t, x + Oi/3Bg) — d 2 x v (t, x)]adj3daB 2 
Jo Jo 
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With the assumption (Hi) we can check that lim^o lE[|-f<s|]^ _1 = 0; from which 
we get d t v(t,x) — G(d% x (t,x)) < 0; hence u is a viscosity supersolution of 
(14. ip . We can analogously prove that u is a viscosity subsolution. But by 
the definition of G-normal distribution 

Bt ~ Vt£ ~ M(®;[g?t,aH}). 
Thus (Bt)t>o is a G-Brownian motion. 

4.1.1 Existence of G-Brownian motion 

In the rest of this course, we denote by Q = C (JS> + ) the space of all M-valued 
continuous paths (uJt)tew.+ with u = 0, equipped with the distance 

OO 

p(w\w 2 ):=y;2- i [(max|a; t 1 -w t 2 |)Al]. 

1=1 

For each fixed T > 0, we consider the following space of random variables: 

(H T ) = L%{F T ) := {X(u) = ^K, • • • ,aO, Vm > 1, 

ii,-,ir»e[0,T] 1 V^C, Lip (r)}. 

It is clear that 7i° C L° p (jF T ), for t <T. We also denote 

oo 
ra=l 



Remark 4.1.4 Obviously Ci.Li P (M m ) and then L^ (Ft), L® p (F) are vector 

lattices. Moreover, since <p, ip G Ciz ji p(W m ) implies <p ■ if) G C;xip(K m ) £/nts 

Ije LJ,CF T ) *mpKe« x • y e ^ P (^r). 

We will consider the canonical space and set B t (u>) = u t , t G [0, oo), for 
u G Q. 

For each fixed T E [0, oo), we set 

L ip (F T ):={<p(B tl ,--- ,B tn ):0<t x ,--- ,t n <T, <p G C LLip (W), n G N}. 

In particular, for each t G [0, oo), B t G L ip (F t ). We are given a sublinear 
function G(a) = Gv,i(a) = \(a + — a 2 a~), a G R. Let £ be a G-normal 
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distributed, or A/"(0; [a 2 , l])-distributed, random variable in a sublinear ex- 
pectation space (fi,?i,E). 

We now introduce a sublinear expectation E denned on H^ = L^JJ^t), 
as well as on H° = L^J 7 ), via the following procedure: For each X G H^ 
with 

X = (p(B tl - B to , B t2 - B tl , ■ ■ ■ , B tm - B tm _ x ) 

for some <p G CiLipiW" 1 ) and = t < t\ < • ■ • < t m < oo, we set 
k[p(B h - B t0 ,B t2 -B tl ,.-., B tm - S tm _J] 



= E[</?(\/£l —t(i€l, • ■ ■ , y/t m — tm-l€m)], 

where (£i, • • • , £ n ) is an m-dimensional G-normal distributed random vector 
in a sublinear expectation space (0,7Y, E) such that £« ~ A/"(0; [a 2 , 1]) and 
such that £ i+ i is independent to (£i, • • • , &) for each i = 1, • • • , m — 1. 

The related conditional expectation of X = <p(B tl , B t , 2 — B tl , ■ • • , B tm — 
Bf m _ 1 ) under "%. is defined by 

E[X|H 4j ] = %(£ tl , fi t2 - B tl , ■ ■ ■ , S tm - S tm _J|W ti ] (4.3) 

= ^(^,■■-,^-5^) 

where 



l/)(xi, • • • , Xj) = E[y?(xi, ■ • • , Xj, ^/tj+l - tjfj+1, • • • , y/t m - *m-lfm)] 

It is easy to check that E[-] consistently defines a sublinear expectation on 
as well as on L® (J 7 ) satisfying (a)-(d) of Definition 13. 1.2L 

Definition 4.1.5 The expectation E[-] : L® (J 7 ) i— > K defined through the 
above procedure is called G— expectation. The corresponding canonical pro- 
cess (B t )t>o in the sublinear expectation space (ft, 7i, E) is called a G-Brownian 
motion. 

Proposition 4.1.6 We list the properties ofK[-\7it] that hold for eachX, Y G 
H* = L%{F): 

(a') IfX > Y, then t[X\H t ] > E[Y\H t ]. 

(b') E[r]\Ht] = V, for each t G [0, oo) and r] G H\. 

(c>) E[X\H t ] - &[Y\H t ] < E[X - Y\H t ). 
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(d>) E[r]X\Ht} = r] + E[X\Ht} + r}-E[-X\H t ], for each r] e H° t . 
We also have 

E [E [X | H t ] \H 8 ]=t[X\ HtAs] , in particular E[E[X\H t ]]=E [X] . 

For each X e L° ip (J*), E[X\H t ] = E[X], where L° ip (J*) = H l T is the linear 
space of random variables of the form 

(p(B t2 - B tl ,B t3 - B t2 , ■ ■ ■ , B tm+1 - B tm ), 
m = 1, 2, • • • , y? e C LLip (R m ), h, ■ ■ ■ , t m , t m+1 , G [t, oo). 

Remark 4.1.7 (b 7 ) and (c') imply: 

E[X + r)\Ht}=E[X\H t }+r)- 
Moreover, ifYe L° ip (F) satisfies E[Y\H t ] = -E[-Y\H t ] then 

E[X + Y\H t ] = E[X\H t } + E[Y\H t }. 
Example 4.1.8 For each s < t, we have E[B t — B S \H S ] = and, for n = 

E[\B t - B s \ n \H s ] = E[\B t _ s \ n ] = . . . / \x\ n exp(-— -)dx. 

y/2ir(t - SJ J-oo l \ l - s ) 

But we have 

E[-\B t - B s \ n \H s \ = E[-\B t _ s \ n \ = -a n E[\B t _ s \ n \. 

Exactly as in classical cases, we have 

E[(B t - B s f\H s \ =t-s, E[(B t - B S ) 4 \H S ] = 3(t - s) 2 , 
E[(B t - B s f\H s \ = 15(t - s) 3 , E[(B t - B S ) 8 \H S ] = 105 (£ - s) 4 , 



n\D, D H \m = w ' 2{t r- s \ n\B t - sjw - 2s/2li " yi ' 2 



TT \/lX 



E[\B t -B s \ 5 \n s ] = 8^ {t " S) ' 2 



71 
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Definition 4.1.9 A process (M t )t>o is called a G '-martingale (respectively, 
G-supermartingale; G-submartingale) if for each t G [0, oo), M t G L® p (Q) 
and for each s G [0, t], we have 

E[M t \H 8 ] = M a , (respectively, < M s ; > M s ). 

Example 4.1.10 (B t ) t > and (—B t ) t > are G-Martingale. (B 2 ) t > is a G- 
submartingale since 

E[B 2 \H S ] = E[(B t - B s f + B 2 S + 2B s (B t - B S )\H S \ 
= E[(B t - B s f] + B 2 = t-s + B 2 > B 2 S . 



40 CHAPTER 4. G-BROWNIAN MOTION: 1-DIMENSIONAL CASE 

4.1.2 Complete spaces of sublinear expectation 

We briefly recall the notion of nonlinear expectations introduced in [55]. 
Following Daniell (see Daniell 1918 [18] ) in his famous Daniell's integration, 
we begin with a vector lattice. 

Let nbea given set and let TC be a vector lattice of real functions defined 
on Q containing 1, namely, TC is a linear space such that 1 G TC and that 
X G TC implies \X\ G TC. TC is a space of random variables. We assume the 
functions on TC are all bounded. Notice that 

aAb = min{a, b} = -(a + b - \a - b\), a V b = -[(-a) A (—&)]. 

Thus X,Y eTC implies that X A Y, X V Y, X + = X V and X~ = (-X)+ 
are all in Ti. 

In this course we are mainly concerned with space H° = L^J 7 ). It 
satisfies 



X u ---,X n en => <p(X ir --,X n )eH, ^eC LL 



rp\ 



Remark 4.1.11 For each fixed p > 1, we observe that 7i^ = {X G Ti, 
E[|X| P ] = 0} is a linear subspace of Ti. To take TCq as our null space, we 
introduce the quotient space TC/TCq. Observe that, for every {X} G TC/Hq 
with a representation X G TC, we can define an expectation K[{X}] := K[X] 

which still satisfies (a) -(d) of Definition \3JJk We set \\X\\ := E[\X\ p }p . 
||-|| forms a Banach norm in TL/TL^. 

Lemma 4.1.12 For r > and 1 < p, q < oo with - + - = 1, we have 

|aH-6| r <max{l,2 r - 1 }(|a| r H-|6| r ), Va,beR; (4.4) 

\a\ p \b\ q 
\ ab \<\^L + \^L. (4.5) 

p q 



Proposition 4.1.13 For each X, Y G L® (J 7 ), we have 

E[\X + Y\ r ] < C r it[\X\ r ] + E[|F| r ]), (4.6) 

E[\XY\] < E[|X| p ] 1/p • E[\Y \ q ] 1/q , (4.7) 

E[\X + Y\ p } 1/p < E[|X| p ] 1/p + E[\Y | p ] 1/p . (4.8) 

In particular, for 1 < p < p' , we have E[|X| p ] 1/p < E[|X| p '] 1/p '. 
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Proof. ( fl~6l) follows from f jOl) . We set 

t X Y 



e[|x|p]Vp e^i?] 1 /? 

By ( 14. 5 p we have 

Eflfnll < E iiL + — < E liL + EPH 
p q p q 

Thus ( 14T7J) follows. We now prove (fl~8l) : 

E[|X + F| p ] = E[|X + Y\ ■ \X + Yf- 1 } 

< E[\X\ ■ \X + F]^ 1 ] + E[|F| • \X + Y^ 1 ] 

< E[|X| p ] 1/p ■ E[|X + Y\ {p - 1)q ] 1/q 
+ E[|F| p ] 1/p • E[\X + Y\ {p - l)q ] l/q . 

We observe that (jp — l)q = p. Thus we have ( 14. 8|) . 
For each p, q > with - + - = 1 we have 

||Xr||=E[|XF|]<||X|| p ||X|| q . 

With this we have ||X|| < ||X|| , if p < p' . ■ 

Remark 4.1.14 It is easy to check thatTi/Ti^ is a normed space under \\-\\ ■ 

We then extend Ti/Ti^ to its completion 7i p under this norm. (7~i p , ||-|L) 

is a Banach space. The sublinear expectation E[-] can be also continuously 
extended from TC/Ho to 7i p , which satisfies (a)-(d). 

For any X 6 TC, the mappings 

X + (lu):H> — >H and I»:Wi — >H 
satisfy 
\X + -Y + \ < \X-Y\ and \X~ - Y~\ = \(-X) + - (-Y) + \ < \X - Y\. 

Thus they are both contraction mappings under ||-|| and can be continuously 

extended to the Banach space (7i p , ||-|L). 

We define the partial order ">" in this Banach space. 
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Definition 4.1.15 An element X in (H, ||-||) is said to be nonnegative, or 
X > 0, < X, if X = X+. We also denote by X >Y, or Y < X, if 
X-Y > 0. 

It is easy to check that X > Y and Y > X implies X = Y in (H p , \\ ■ || ). 

The sublinear expectation E[-] can be continuously extended to (H p , ||-|L) 
on which (a)— (d) still hold. 

4.1.3 G-Brownian motion in a complete sublinear ex- 
pectation space 

We can check that, for each p > and for each X G L° t JT) with the form 
X{u) = <p(B tl ,- ■ ■ , BtJ, for some <p G C LLip (R m ), 

E[|X|] = ^^ E[|X| p ]=0 ^^ ^)E0,Vxel m 

For each p > 1, ||X|| p := E[|X| P ]?, X G L° ip (F T ) (respectively, L%(F)) 
forms a norm and that L® (J 7 ?) (respectively, L^J 7 )), can be continuously 
extended to a Banach space, denoted by 

Ht = L p g (!Ft) (respectively, TC = LP G (T)). 

For each 0<£<T<oowe have L v G {Tt) C L g {Tt) C L v g {T}. It is easy 
to check that, in L g {Tt) (respectively, Lq(.F)), E[-] still satisfies (a)-(d) in 
Definition [3X1 

We now consider the conditional expectation introduced in (j4.3p . For 
each fixed t = tj < T, the conditional expectation E[-|7i t ] : L° p (jF T ) i— »■ 7i° is 

a continuous mapping under ||-|| since E[E[X|7i t ]] = E[X], X G L° p (jF T ) and 
E[E[X|ft t ] - E[K|7it]] < E[X - Y], 



E[X\H t ] - E[Y\H t 



< \\X -Y\ 



It follows that E[- \Ttt] can be also extended as a continuous mapping L p g (Tt) 
L v G [Tt). If the above T is not fixed, then we can obtain E[-|7i t ] : L G (T) y 



Proposition 4.1.16 The properties of Proposition 4-1-6 ofE[-\7i t ] still hold 
forX,YeL G (F): 
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(a') IfX > Y, then E[X\H t ] > E[Y\H t }. 

(b') E[r)\Ht\ = V, for each t E [0, oo) and r) E L G (T^. 

(c') E[X\H t ] - E[Y\H t ] < E[X - Y\H t }. 

(d>) E{rjX\H t ] = ri + E[X\H t } + r]-E[-X\H t ], for each bounded q E L l G (F t ). 

(e>) E[E[X\H t ]\H s } = E[X\H tAs }, in particular, E[E[X\H t ]] = E[X]. 

(f>) For each X E L X G (J%) we have E[X\H t ] =E[X]. 

Definition 4.1.17 An X E L G (J-) is said to be independent of Tt under 
the G -expectation E for some given t E [0, oo), if for each real function $ 
suitably defined on M. such that $(X) E L^(jF) we have 

E[$(X)|ft t ]=E[$(X)]. 



Remark 4.1.18 It is clear that all elements in L^F) are independent of 
JF . Just like the classical situation, the increments of G-Brownian motion 
(B t+S — B s ) t > is independent of T s . In fact it is a new G-Brownian motion 
since, just like the classical situation, the increments of B are identically 
distributed. 

Example 4.1.19 For each n E N, < £ < oo and X E L G {!Ft), since 
E[B^S t 1 } = El-B^Tt 1 ], we have, by (f) of Proposition ^ 1.1 6\ 



E[X(B T - B t ) 2n - 1 } = E[X + E[(B T - B t ) 2n - l \H t ] + X'E[-(B T - B t ) 2n ^ 1 \'H^ 

= E[\X\]-E[B 2 T n - 1 ], 
E[X(B T - B t )\H t \ = E[-X(B T - B t )\H t ] = 0. 

We also have 

E[X(B T - B t ) 2 \H t } = X + (T -t)- a 2 X-(T - t). 



Remark 4.1.20 It is clear that we can define an expectation E[-] on L® (J 7 ) 
in the same way as in Definition \4-l-l\ with the standard normal distribution 
F =JV(0, 1) m place o/F 5 = JV(0; [a, 1]) on (R, B(R)). Since F is dominated 
by F^ in the sense F[(p] — F [ip] < ¥^[cp — ip\, then E[-] can be continuously 
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extended to L G (J-). E[-] is a linear expectation under which (B t )t>o behaves 
as a Brownian motion. We have 

E[X] < E[X], VX e L G {F). (4.9) 

In particular, E[Bl n S t 1 } = Ef-5^ 1 ] > Ef-5^ 1 ] = 0. Such kind of exten- 
sion under a domination relation was discussed in details in [55]. 

The following property is very useful. 

Proposition 4.1.21 LetX,Y e L G (T) be such that E[Y\H t } = -E[-Y\H t ], 
then we have 

E[X + Y\H t ] = E[X\H t ] + W\H t }. 

In particular, if t = and E[Y\H ] =t[Y] = E[-Y] = 0, then E[X + Y] = 

E[X]. 

Proof. We just need to use twice the sub-additivity of E[-|7i t ]: 

E[X + Y\H t ] > E[X\H t ] - E[-Y\H t ] = E[X\H t ] + E[Y\H t ] 
>E[X + Y\H t ]. 



Example 4.1.22 We have 

E[B 2 - B 2 S \H S ] = E[(B t -B s + B s ) 2 - B 2 S \H S ] 

= E[{B t - B s ) 2 + 2{B t - B S )B S \H S 
= t — s, 



since 2(B t — B S )B S satisfies the condition for Y in Proposition \4-l- 21\ and 



E[{B 2 t - B 2 s y\n s ] = E[{(B t -b s + B s y - B 2 s y\n s ] 

= E[{(B t - B s f + 2(B t - B S )B S } 2 \H S ] 
= E[(B t - B s ) 4 + A(B t - B s fB s + A(B t - B S ) 2 B 2 \H S ] 
< E[(B t - B s ) 4 ] + AE[\B t - B S \ 3 ]\B S \ + 4(t - s)B 2 
= 3(£ - s) 2 + 8(t - s) 3/2 \B s \ + 4(t - s)B 2 . 



4.1. 1-DIMENSIONAL G-BROWNIAN MOTION 45 

4.1.4 Ito's integral of G— Brownian motion 

Bochner's integral 

Definition 4.1.23 For T e M. + , a partition ttt of [0,T] is a finite ordered 
subset Ti = {t%, ■ ■ ■ , t^} such that — to < ti < • ■ ■ < tN — T . We denote 

/i(7r T ) = max{|t i+ i -U\: i = 0,1,- ■ ■ ,N - 1}. 

We use 71? = {to,ti,--- ,t^} to denote a sequence of partitions of [0,T] 
such that linijv^oo ^(^t ) = 0- 

Let p > 1 be fixed. We consider the following type of simple processes: 
For a given partition {to, ■ • • , t^} = tit of [0, T], we set 

N-l 

m{w) = 5Z^Mi[t J; * i+1 )(*), 

where & G L G {Tt^), i = 0, 1, 2, ■ • ■ , N — 1, are given. The collection of this 
form of processes is denoted by M G (0, T). 

Definition 4.1.24 For an 77 e Mg'°(0,T) with Vt = En^wlIfe^oW 
i/ie related Bochner integral is 

N-l 



■/o i=0 



P,0/ 



Remark 4.1.25 For eac/j 77 e Mg u (0,T) we set 

rp N—l 

E T [ V ] := i / %]d* = - J] E[^H](t i+1 - £,-)• 
•>o i=0 

It is easy to check that E^ : M§ (0,T) 1 — > R forms a sublinear expectation 
satisfying (a) -(d) of Definition ^. l.H From Remark 4. 1.11 we can introduce 



a natural norm \\v\\m p (ot) = 1 Jo ^ilVtl^dt \ . Under this norm M G (0, T) 
can 6e continuously extended to a Banach space. 
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Definition 4.1.26 For each p > 1, we will denote by M G (0, T) the comple- 
tion of Mq (0, T) under the norm 

IMUg(o,T) = \J &[\Vt\ p ]dt> . 
For T) G Mq (0, T), we have 

\^ j\ t (co)dtf] < ^nu^nt 1+ i - h 

E[\Vt\ p ]dt. 



, fT n N-l 

E[ 

1 ' T 



r./ 



We then have 



Proposition 4.1.27 27ie linear mapping J r]t(uj)dt : Mq (0, T) i— > L g {Tt) 
is continuous; thus it can be continuously extended to M G (0,T) i— >■ L g [Tt)- 
We still denote this extended mapping by j Q i] t (uj)dt, rj G M G (0, T). VFe /iave 



cT 1 /-T 



E[|^y Vt(w)dm<±J n\Vt\ p ]dt, \/ V eM G (0,T). (4.10) 

We have M G (0, T) D M G (0, T), for p < q. 

Ito's integral of G— Brownian motion 

Definition 4.1.28 For each r] G M G '°(0,T) with the form 

N-l 
3=0 



we define 



T N-l 

I{ V )= / V (s)dB s :=J2^(B tj+1 -B tj ). 
Jo j=0 
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Lemma 4.1.29 The mapping I : M G (0, T) i — > L g [Tt) zs a linear con- 
tinuous mapping and thus can be continuously extended to I : M G (0,T) i — > 
1? g {Tt)- We have 

E[[ r l (s)dB s ]=0, (4.11) 

E[( / V (s)dB s ) 2 ] < f E[( V (t)) 2 ]dt. (4.12) 
Jo Jo 



Definition 4.1.30 We define, for a fixed T] G M g (0,T), the stochastic inte- 
gral 

I r){s)dB s ;=I{ri). 

Jo 

It is clear that $U$ and (g[7|) still hold for rj G M G (0, T). 



Proof of Lemma 14.1.291 From Example 14.1.191 f° r eacri 3, 
E[^(B tj+1 - B t .)\H tj ] = E[-^(B t . +1 - B t .)\H tj ] = 0. 
We have 

E[ / V (s)dB s ] = E[ / r]{s)dB s + ^ x {B tN - B tN ^)} 
Jo Jo 

ptN-1 

= E[ / r]{s)dB s + E[£ N „ 1 {B tN - Bt^JlHt^]] 
Jo 

/•tjv-i 

= E[ / r]{s)dB s ]. 
Jo 

We then can repeat this procedure to obtain (14.111) . We now prove (14.121) : 
E[(jf ri(s)dB a ) 2 ]=&[n"~\(s)dB a + Z N _ 1 {B t „-B tN _S\ ] 

/ 77(s)rfS s J +E[2( / »7( S )dSj ^-i(Bt w - St^; 



+ ^N-l( B t N ~ B tN _ 1 ) |?4v-l]] 

= E[(/ \(s)dB s ) +t; 2 N _ 1 (t N -t N - 1 )]. 



o 
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Thus E[(/ tjv v(s)dB s ) 2 ] < E[ ( J^" 1 77(5)^) *] + E&^J^ - ^_i). We 
then repeat this procedure to deduce 



T TV— 1 t 

E[(/ i7(s)dB.) 2 ]<X;^) 2 ](*i + i-*i)=/ E[(^(t))Vt. 

>>0 „_ n ./0 



We list some main properties of the Ito's integral of G-Brownian motion. 
We denote for some < s < t < T, 



/ r] u dB u := / I[ s , t ](u)r) u dB u . 



We have 



Proposition 4.1.31 Let rj, 9 6 Af£(0,T) and let < s < r < t < T. Then 

in Lq(J-'t) we have 

(i) f s r] u dB u = f s r) u dB u + f r rj u dB u , 

(ii) j s (ai] u + 6 u )dB u = a j s i] u dB u + j s 6 u dB u , if a is bounded and in L^F), 

(Hi) E[X + f^rj u dB u \H s } = &[X], VX e L l G (F). 



Quadratic variation process of G— Brownian motion 

We now study a very interesting process of the G-Brownian motion. Let n 
N — 1, 2, • • • , be a sequence of partitions of [0, t]. We consider 



N 

t ' 



AT-l 



B t = J2^ B t? +1 - B t? 



f tN \ 

"3 + 1 l 3 

3=0 

AT-l JV-1 

= Yl 2B t?( B t? +l - B tf) + X^f+i ~ B tf ) 2 - 
3=0 3=0 

As fJ.{7r^) — > the first term of the right side tends to 2 J B s dB s . The second 
term must converge. We denote its limit by (B) t , i.e., 



N-i r t 

(B) f = lim y^(B t N -B,n) 2 = B?-2 B s dB s . (4.13) 
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By the above construction {{B) t } t > , is an increasing process with (B) = 0. 
We call it the quadratic variation process of the G-Brownian motion B. 
Clearly (B) is an increasing process. It characterizes the part of statistic 
uncertainty of G-Brownian motion. It is important to keep in mind that 
(B) t is not a deterministic process unless the case a = 1, i.e., when B is a 
classical Brownian motion. In fact we have 

Lemma 4.1.32 We have, for each < s < t < oo 

E[(B) t -{B),\7i.]=t-s, (4.14) 

E[-((B) t -(B) s )\n s ] = -a\t-s). (4.15) 



Proof. By the definition of (B) and Proposition I4.1.3"l] -(iii). 
E[(B) t - (B) s \H S ] = E[B 2 - B 2 - 2 f B u dB u \H s ] 

J s 

= E[B? - B 2 S \H S ] =t-s. 

The last step can be checked as in Example 14.1.221 We then have (14.141) . 
( 14.151) can be proved analogously with the consideration of E[— (B 2 — B 2 ) \H S ] = 
-a 2 (t-s). m 

To define the integration of a process rj 6 Mq(0, T) with respect to d (B), 
we first define a mapping: 

„ T N-l 

QoAv) = / v(s)d(B) s := J>«B) t . +1 - (B) t .) : M^°(0,T) -> L\T T ). 

Lemma 4.1.33 For each r] e M^'°(0,T) ; 

n\QoAv)\] < [ n\Vs\}ds. (4.16) 

Jo 

Thus Q ,t : Mq (0,T) i— > L 1 (JF T ) zs a continuous linear mapping. Con- 
sequently, Q ,t can fre uniquely extended to L^(0,T). We still denote this 
mapping by 

r r ] (s)d(B) s = Qo, T (v), r/eM^(0,T). 



/o 
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We still have 

E[| / rj(s)d(B) s \}< [ E[\ Vs \]ds, V V eM G (0,T). (4.17) 

Jo Jo 

Proof. From Lemma [4.1.321 (14.161) can be checked as follows: 

AT-l N-l 

j=0 
JV-l 

j=0 
T 

t[\n,\]ds. 



A very interesting point of the quadratic variation process (B) is, just like 
the G-Brownian motion B itself, the increment (B) t+s — (B) s is independent 
of T s and identically distributed like (B) t . In fact we have 

Lemma 4.1.34 For each fixed s > 0, ((B) s+t — (B) s ) t > is independent 
of T s . It is the quadratic variation process of the Brownian motion B^ = 
B s+t -B s ,t> 0, i.e., (B) s+t - (B) s = (B% 

Proof. The independence follows directly from 

(B) s+t - (B) s = B 2 t+S - 2 f B r dB r - [Bl - 2 f B r dB r ] 

Jo Jo 

/s+t 
(B r — B s )d(B r — B s ) 

■ 

Proposition 4.1.35 Let < s < t, £ G L l G (F s ). Then 

E[X + i(B 2 t - Bl)} = E[X + £(B t - B s ) 2 } 

= E[X + a(B) t -(B) s )]. 
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Proof. By ( 14.131) and Proposition 14.1.211 we have 

E[X + H{B 2 t - B 2 )] = E[X + a(B) t - (B) s + 2 f B u dB u )} 

J s 

= E[X + a(B) t -(B) s )]. 
We also have 

E[X + £(B 2 - B 2 S )} = E[X + £{(Bt - B s f + 2(B t - B S )B S }] 

= E[X + £(B t -B s ) 2 }. 

■ 

We have the following isometry 

Proposition 4.1.36 Let 7] e M G (0,T). We have 

E[( !\{s)dB s f] = E[ frf{s)d (B) s \. (4.18) 

jo Jo 



Proof. We first consider rj G Mq (0,T) with the form 

N-l 

m(u) = 5^6-( w ) I fe,*i+i)(*) 

i=o 

and thus f Q r](s)dB s := J2j=o £,j{Bt i+1 — B^). From Proposition 14.1.211 we 
have 

E[X + 2^(B tj+1 - B t] )UB u+1 - B u )\ = E[X], for X e L G (F), i + J- 

Thus 

r T /N-l \ 2 AT-1 

E[(/ v(s)dB s ) 2 }=E[[J2UBt J+1 -B tj )) ]=E[J2%(Bt j+1 -B tj f]. 

From this and Proposition 14.1.351 it follows that 

E[(/ v(s)dB s ) 2 ]=E[J2^((B) t]+1 - (B) tj )]=n f v\s)d(B) s ]. 
Jo j=0 Jo 

Thus ( I4.18P holds for 77 e M G ' (0, T). We can continuously extend the above 
equality to the case r\ e M^(0,T) and prove (j4.18p . ■ 
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The distribution of (B) t 

The quadratic variation process (B) of G-Brownian motion B is a very in- 
teresting process. We have seen that the G-Brownian motion B is a typical 
process with variance uncertainty but without mean-uncertainty This uncer- 
tainty is concentrated in (B). Moreover, (B) itself is a typical process with 
mean- variance. This fact will be applied to measure the mean-uncertainty of 
risky positions. 

Lemma 4.1.37 We have 

E[{B s ) 2 t \H S ] = E[(5) 2 ] < 10t 2 . (4.19) 



Proof. 



E[(B) 



2 }.=E[{{B t f- 


2 / 

Jo 


B u dB t 


y] 


< 2E[(B t ) 4 } + 8E[( 


/ B u dB u ) 
Jo 


< 6t 2 + 8 / 
Jo 


ip., 


■j) 2 ]du 




= 10t 2 . 









Definition 4.1.38 A random variable £ of a sublinear expectation space 
(f2,7i, E) is said to be 1A\p,M- distributed for some given interval [fJ,,~p] C R 
if for each ip G CYlj P (R) we have 

Efa(0] = SU P_ v( x )- 
xe[/i,/x] 



Theorem 4.1.39 Let (b t ) t > be a process defined in (Q,7i,E) such that 

(i) b = 0. 

(ii) For each t, s > 0, the difference b t+s — b t and b s are identically distributed 

and independent to (b tl , b t2 , • ■ • , b tn ) for each n GN and < t±, • • ■ , t n < t. 

(Hi) b = and lim no i[&f]t -1 = 0. 

Then b t is U^t^t}- distributed with /I=E[6i] and \x = — E[— bi]. 



4.1. 1-DIMENSIONAL G-BROWNIAN MOTION 53 

Proof. We need only to prove that b t is ^^-distributed. We first prove 

that 

E[b t ] = Jit and - E[-b t ] = fxt. 

We set ip(t) := E[b t ]. Then (p(0) = and lim^ (fit) — > 0- Since for each 
t,s>0 

cp(t + s) = E[bt +S ] = E[(b t+S - b s ) + b s ] 
= (p(t)+(p(s). 

Thus (f(t) is linear and uniformly continuous in t which means E[b t ] = Jit. 
Similarly — E[— b t ] = lit. 

We now prove that b t is U^-pt] -distributed. We just need to prove that 
for each fixed tp e CiL ip {R), the function 

u(t,x) :=E[ip(x + b t )}, (t,x) E [0,oo) xl 
is the viscosity solution of the following G Mi p-drift equation 

d t u - 2G{d x u) = 0, for t > 0, u\ t=0 = if (4.20) 

with G(a) = Gfj, t j[(a) = \(jla + — /ia _ ). Since the function 

u(t, x) = max if(x + vi) 

is the unique solution of the PDE (14.201) . 

We first prove that u is locally Lipschitz in (t, x). In fact, for each fixed 
t, u(t, •) eC LLip (R) since 

\E[ip{x + bt)] - E[<p{y + b t )]\ < \E[\ip{x + h) - <p{y + h)\]\ 

<\E[C(l + \x\ m + \y\ m +\bt\ m )\x-y\]\ 
< C t (l + \x\ m + \y\ m )\x - y\. 

For each 5 G [0, t], since bt — bs is independent to bs, we also have 

u(t,x)=E[(p(x + bs+(b t -b 5 )] 

= E[E[ip(y+(b t -b s ))] y=x+bs } 
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hence 

u(t,x)=E[u(t-5,x + b 5 )]. (4.21) 

Thus 

\u(t, x) — u(t — S,x)\ = \E[u(t — 5,x + b$) — u(t — 5, x)] \ 

< E[\u(t-S,x + b s ) -u(t-5,x)\] 

<E[C(l + \x\ m +\b 5 \ m )\bs\] 

< d(l + \x\ m )8. 

To prove that u is a viscosity solution of the PDE (14.201) we fix a time- 
space point (t,x) G (0, oo) x M. and let v G C b ' ([0, oo) x E) be such that 
v > u and v(t,x) = u(t,x). From ( 14.211) we have 

v(t, x) = E[u(t -5 } x + b s )\ < E[v(t -5,x + b s )\ 

Therefore from Taylor's expansion 

< E[v{t -S,x + b§)- v(t, x)] 
= E[v(t -S,x + b s )- v(t, x + b s ) + (v(t, x + b 5 )- v(t, x))} 
= E[-d t v(t, x)5 + d x v(t, x)b s + Is] 
< -d t v{t, x)5 + E[d x v(t, x)b s ] + E[I S ] 
= -d t v(t, x)S + 2G(d x (t, x))S + E[I 5 }. 

where 

h= [ l-d t v(t-f35,x + (3bs) + d t v(t,x)}dp5 
Jo 

+ / [d x v(t - p5,x + f3b s ) - d x v(t,x)]d(3b s . 
Jo 

With the assumption that limtjoEf&^t -1 = we can check that 

limE[|/ (5 |]r 1 = 0; 

<5|0 

from which we get d t v(t, x) —2G(d x (t, x)) < 0; hence u is a viscosity superso- 
lution of (I4.20p . We can analogously prove that wis a viscosity subsolution. 
It follows that b\ is U^-m -distributed. The proof is complete. ■ 
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Corollary 4.1.40 For each t <T < oo ; we have 

fMT-t) <b T -b t <Jl(T-t), mL^(T). 
Proof. It is a direct consequence of 

E[(b T -bt-(T- £))+] = sup_E[(r? - Jl) + (T - t)] = 

and 

E[(b T -h- o~ 2 (T - t))~] = sup_E[(r] - fjT)-(T - t)] = 0. 

/J,<ri<Ji 

■ 

Corollary 4.1.41 We have 

E[(B S ) 2 \H S ] = E[(B) 2 ] = t 2 (4.22) 

as well as 

E[(B°f t \H S ] = E[(B} 2 ] = t\ E[(B*)t \H S ] = E[(B) 4 t ] = t\ 

Theorem 4.1.42 For each x G R, Z G M%(Q,T) and r] G M^(0,T) the 
process 

M t = x+ [ Z s dB t + I r] s d (B) s - [ 2G(r) s )ds, t G [0,T] 
Jo Jo Jo 

is a martingale: 

E[M t \H s ] = M s , 0<s<t<T. 

Proof. Since E[J s Z r dB t \H s ] = E[— J s Z r dB t \H s ] = we only need to prove 

that 

?t ft 



M t = [ Vsd(B} s - [ 2G( Vs )ds 
Jo Jo 



is a G-martingale. It suffices to consider the case where 77 is a simple process, 
i.e., rjt = ^2 k = £fcl[*k,tfc+i)(0- m f ac t we only need to consider one-step case 
in which we have 

E[M tk+1 - M tk \H tk ] = EM(B) tk+1 - (B) tk ) - 2G(^ k )(t k+1 - t k )\H tk ] 
= EM(B) tk+1 - (B) tk )\H tk ] - 2G(H k )(t k+1 - t k ) 
= Cn(B} tk+1 ~ (B) tk ]+^E[-((B) tk+i - (B) tk )\ 

— 2G(^ k )(t k+ i — t k ) 
= 0. 
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Problem 4.1.43 For each £ = (f{B tl — B to , • • • , B tN — B tN _ l ) we have the 
following representation: 

e = E[£]+ f z s dB s + [ Vsd(B) s - [ 2G(rj s )ds. 
Jo Jo Jo 

4.2 Ito's formula for G— Brownian motion 

We have the corresponding Ito's formula of $(X t ) for a "G-Ito process" X. 
For simplification we only treat the case where the function $ is sufficiently 
regular. We first consider a simple situation. 

Lemma 4.2.1 Let $ G C 2 (W l ) be bounded with bounded derivatives and 
{d^ X "^}^,u=i are uniformly Lipschitz. Let s G [0, T] be fixed and let X = 
{X 1 , ■ ■ ■ , X n ) T be an n-dimensional process on [s, T] of the form 

XI = X: + oT{t -s) + V "((B) t - (B) s ) + f3»(B t - B s ), 

where, for v — 1, • • • ,n, a u , rf and (3 U , are bounded elements o/L^,(jF s ) and 
X s = (Xj, • • • ,X™) T is a given M. n -vector in L^.(jF s ). Then we have 



$(X t ) - $(X S ) = / d x u$(X u )P"dB u + d x MX u )a v du (4.23) 

J s J s 

+ J\d x ^(X u ) V » + l -dl» x MXu)^nd (B) u . 

Here we use the Einstein convention, i.e., each single term with repeated 
indices /i and/or v implies the summation. 

Proof. For each positive integer N we set 5 = (t — s)/N and take the 

partition 



7T, 



s,t] 



{t«,t? r --,t N N } = {s,s + 5 r --,s + N5 = t}. 



We have 

N-l 



k=0 

N-l 



$(X S ) + J2[d x »<f>(X t N)(X? N - xjt 

■*— ' k t fe+l l k 



k=0 

1 



+ ^[dl^{X tN ){X» N - Xfr){X% - X£) + rtf]] (4.24) 

Z k fc+l k fc+1 k 
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where 

Vk = [a&wHXtg +e k (x t ^ i -x^ ))-d 2 x , x ^(x t ^)}(x^ + -x^)(x^ +r x^) 

with O k E [0,1]. We have 



E[|r/f |] = E[\[dl x „<f>(x t » + e k {x q - x tN )) - d^ x Mx t »)](x? N ~ X?* 

k fc+1 k k T- k+1 *> k 



x 



<cE[|X Ci --X 4f | 3 ]<C[5 3 + 5 3 / 2 ], 



where c is the Lipschitz constant of {d xlix u^}^ u=1 - Thus ^ fc E[|?7^|] — > 0. 
The rest terms in the summation of the right side of (14.241) are £^ + (^ with 



JV-l 



+ l -dl» x MX t *WP{B tNk+i - B t M)(B t » +i - Bq)} 
and 

JV-l 



cf = \ E ^*(%)K(Ci - tf ) + *A<^ +1 - (5)^)1 

fc=0 

xK(Ci-^) + ^((%-(5) tf )] 

fc + 1 A: 

+ ^K(if +1 - tf ) + 7A(s) tf+i - (s) tf )](s tf+i - B q ). 

We observe that, for each u G [i^, ^fc+i) 



JV-l 



E[|^$(X U ) - X^*(Xtf)I [t *Wu)| 2 ] 



fc=0 



<^i 2 ] 



k 

<c 2 E[\X u -X tN \ 2 } <C[5 + 5 2 ] 

h 
^JV-l 



Thus ^^T i ^$(X^)I [ ^ i ^ +i) (-) tends to <9^$(X) in M^(0,T). Similarly, 

JV-l 

Y,%» x MX t »)\t»,t»){-) -* dl x „$(X.) in M 2 (0,T). 



fc=0 
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Let N — > oo. From the definitions of the integrations with respect to dt, dB t 
and d (B) t the limit of £f in L%{T^) is just the right hand of (jQ5l> . By the 
estimates of the next remark we also have ^ — > in L^.(jF 4 ). We then have 
proved (1423]) . ■ 



Remark 4.2.2 VFe nai>e the following estimates: For ip e M G ' (0, T) sitcn 
£na£ Vf = EhC 1 !?.^!)^) anrf ^t = {° < t 0l --- ,t N = T} with 
lim^oo MO = and ^f^ 1 E[|^|](tf +1 - tf) < C /or aU N = 1, 2, • • • , 

we nave 

N-l 



E[ix;^(d-o 2 i]-o, 



fc=0 



and, thanks to Lemma 4-1-34 



N-l N-l 



fc+i 
fc=0 A:=0 



k=0 

N-l 

E^iefiKCi-^ 



as we// as 



JV-l 



N-l 

< Y,mk\n((B) t M +i - (B) t »)\B t » +i - B t ,\\ 

fc=0 
N-l 

< T.^klMiB)^ - (B) t ^E[\B t , +i -B t ,\^ 

fc=0 

N-l 

= £%?|](&i-tf) 8/2 ->o. 



fc=0 
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We also have 

N-l 

E[]E^((^ +1 -(^)(Ci-^)l] 

fc=0 
N-l 

< £%f KCi -tk) -n((B) t , +i - (B) t ,)\n t ,]} 

k=0 

N-l 



fe=0 



and 

N-l N-l 



fc=0 k=0 



= v^ i[| ^ |](t ^-^ )3/2 ^°- 

■ 
We now consider a more general form of Ito's formula. Consider 

X? = XZ+ f a v 8 ds + [ V :d(B,B) s + [ (3»dB s . 
Jo Jo Jo 

Proposition 4.2.3 Let a v , [3 U and r\ v , v — 1, • • • ,n, are bounded processes 
of Mq(0, T). Then for each t > and $ in L 2 G {Tt) we have 

$(X t ) - $(X 8 ) = I d x »$(X u )(3»dB u + / d x MX u )a u u du (4.25) 

J s J s 



[d x ^(x u ) V : + l -di, x MXu)^:]d (b) u . 



Proof. We first consider the case where a, i] and f3 are step processes of 
the form 

N-l 

•rhiu) = ^2tk{u)l[t k ,t k+1 ){t). 

fc=0 
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From Lemma [4.2.11 it is clear that ( 14. 25ft holds true. Now let 

X? N = X» + fctfda+ f\:' N d(B) s + f ^ N dB s , 
Jo Jo Jo 

where a N , r] N and f3 N are uniformly bounded step processes that converge 
to a, 7] and j3 in Mq(0, T) a.s N —>■ oo, respectively. From Lemma [4.2.11 

$(X? N ) - $(X ) = f * d a r*(X?)!% lf dB u + f d x MXu)<l N du (4.26) 



+ j\d x MXu)rt N + \dl» x MXuWu N P V u N \d (B) u . 

Since 

E[\X? N - X»\*\ < 3E[| /"(of - a s )ds\ 2 ] + 3E[| f (^ N - rftd (B) s | 2 ] 

Jo Jo 

+3E[| [\(3:> N -f3:)dB s \ 2 ]<3 [ &[(<#* - off]ds + 3 / E[|^ - tf| 2 ]<fc 
Jo Jo Jo 

+3 f tm N - (3:y\ds. 

Jo 



Furthermore 



d x ^(X N )r] v ' N + d 2 x » x MX N )^ N P' N -> ^$(X)^ + ^^(X)/^/?' 

^(X.^K'* - ^$(X.K; 
d^X^p?" -* d x »$(X.)(3 v . 

We then can pass limit in both sides of (14.261) and get (03 



4.3. STOCHASTIC DIFFERENTIAL EQUATIONS 61 



4.3 Stochastic differential equations 

We consider the following SDE defined on M&(0, T; R n ): 

X t = X + [ b(X s )ds+ [ h(X s )d(B) s + [ a(X s )dB s , te [0,T]. (4.27) 
Jo Jo Jo 

where the initial condition Xo G M. n is given and b,h,a : M n h-> R n are given 

Lipschitz functions, i.e., |y(a;) — <£>(a/)| < iiT|x — s'|, for each x, x' G M n , 

if = b, h and a, respectively. Here the horizon [0, T] can be arbitrarily large. 

The solution is a process X G M^(0,T;M n ) satisfying the above SDE. We 

first introduce the following mapping on a fixed interval [0, T]: 

A.(Y) := F G M^(0,T;M n ) i — ► M^(0,T;R") 

by setting A t with 

A t (Y) =X + [ b(Y s )ds + [ h(Y s )d(B) s + [ a(Y s )dB s , £G [0,T]. 
Jo Jo Jo 

We immediately have 

Lemma 4.3.1 For each Y,Y' G Mg(0,T;R") we have the following esti- 
mate: 

E[|A t (y)-A t (r)| 2 ]<C ftHY-Y^ds, tG[0,T], 

Jo 

where C = 3K 2 . 

Proof. This is a direct consequence of the inequalities (14.101) . (14.121) and 

(EH). ■ 

We now prove that SDE (I4.27P has a unique solution. By multiplying 
e -2Ct on both sides Q f the above inequality and then integrating them on 
[0, T] it follows that 

/ E[|A t (r) - A t (Y')\ 2 }e- 2Ct dt < C [ e- 2Ct [ E[\Y S - Y^\ 2 }dsdt 
Jo Jo Jo 

i-T r T 



C f f e~ 2Ct dtE[\Y s - Y^\ 2 ]ds 

(2C)~ 1 C / V 2Cs - e- 2CT )E[\Y s - Y' s \ 2 ]ds. 
Jo 
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We then have 

/ T E[|A t (F) - A t (Y')\ 2 ]e- 2Ct dt < ~ f ' E[\Y t - Y*\ 2 ]e-**dt. 
Jo ^ Jo 

We observe that the following two norms are equivalent in Mq(0, T; M n ): 

/ E[\Y t \ 2 ]dt~ f E[\Y t \ 2 ]e- 2Ct dt. 
Jo Jo 

From this estimate we can obtain that A(Y) is a contraction mapping. Con- 
sequently, we have 



Theorem 4.3.2 There exists a unique solution X G Mq(0, T;R n ) of the 
stochastic differential equation ( f^.£7| ). 



4.4 Backward SDE 

We consider the following type of BSDE 

Y t = E[^ + J f(s,Y s )ds\H t }, te[0,T}. (4.28) 

where £ G L^(T T ; M m ), /(t, y) G M^(0, T; M m ), t e [0, T], y G M"\ are given 
such that 

|/(t,y)-/(t,i/0|<*|y-i/|. 
We have 

Theorem 4.4.1 There exists a unique solution (it)te[o,Tl £ M 1 (0, T). 
A t (Y)=E[£+ / /(s,y s )ds|W t ] 
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f T E[\A t (Y) - A t (Y')\]e^dt < f t[ f \f(s,Y s ) - f(sX)\ds}e^dt 
o Jo Jt 

<C / / E[|y a -y a '|]e^dsdt 
Jo A 

= C / E[|Y,-y,'|] / e^dtds 
Jo Jo 

= jJ o T n\Y s -YMe^-l)ds 



<£ /" E[|y, -Yl'lle^ds. 



■s - 1 s 



0./O 

W^e choose f3 = 2C, then 

f T E[\A t (Y) - A t (Y')\]e^dt < \ f t[\Ys ~ Y^ s ds. 
Jo * Jo 

Excercise. We define a deterministic function u = u(t, x) on (t, x) G [0, T] x 

R d by 

dXf = f(Xl' x )ds + h(Xl' x )d (B) s + a(Xl' x )dB s , s e [t, T], Xf = x 

where $ and (7 are given IR m - valued Lipschitz functions on IR d and IR d x R m , 
respectively, w is a viscosity solution of a fully nonlinear parabolic PDE. Try 
to find this PDE and prove this interpretation. 
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Chapter 5 

Vector Valued G-Brownian 
Motion 

5.1 Multidimensional G— normal distributions 

For a given positive integer n we will denote by (x, y) the scalar product of 
x, y G W 1 and by \x\ = (x, x) 1 ^ 2 the Euclidean norm of x. We denote by 
§(n) the space of all d x d symmetric matrices. §(n) is a 2 -dimensional 
Euclidean space of the norm (A,B) :=tr[AB]. We also denote E>(n) + subset 
of non- negative matrices in S(n). 

For a given rf-dimensional random vector £ = (^, • • • , £<$) in a sublinear 
expectation space (f2, 7i, E) with zero mean uncertainty: 

Efc] = (E[£i],--.,E[&]) = 0, EK] = 0. 

We set 

G(A) = Gt(A):=±&[(M,t)]> A em- 

This function G(-) : §(n) i — ► R characterizes the variance-uncertainty of £. 
Since £ G 7i implies E[|£| 2 ] < oo, there exists a constant C such that 

|G(A)| <C\A\, VAG%). 



Proposition 5.1.1 ITie function G(-) : §(n) i — > M is a monotonia and 
sublinear mapping: 

65 
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(a) A>B =► G{A) > G(B); 

(b) G(XA) = \ + G(A) + \-G(-A); 
(c)G{A + B) < G{A) + G{B). 

Moreover, there exists a bounded, convex and closed subset V C E> + (d) (the 
non-negative elements ofS(d)) such that 

G(A) = \ sup (A, B). (5.1) 

* Ber 

We introduce rf-dimensional G-normal distribution. 

Definition 5.1.2 In a sublinear expectation space (0,E, TC) a d- dimensional 
random vector £ = (£i, • • • ,^) e 7i d with is said to be G-normal distributed 
if for any random vector ( in 7i d independent to £ such that ( ~ £ we have 



< + b( ~ Va 2 + 6 2 £, Va, 6 > 0. 
i/ere G(-) zs a real sublinear function defined on S(n) by 

G(A):=±t[{M,t)], AeS{n). 

Since G is uniquely determined byTc. S(n) via Ii5.1\) £ zs a/so called jV(0, T) 
distributed, denoted by £ ~ A/"(0, T). 

Proposition 13.2.71 tells us how to construct the above (£, C) based on the 
distribution of £. 

Remark 5.1.3 By the above definition, we have v^Ef^] = E[£ + C] = 2E[£] 
and \/2E[-£] = E[-f - (] = 2E[-£] it /o//ows t/iat 

E[e] = i[-e] = o, 

Namely G-normal distributed random variable £ has no mean uncertainty. 

Remark 5.1.4 If C, is independent to ( and £ ~ (, such that A3.2fy satisfies, 
then — £ is independent to —( and — £ ~ —(. We also have a(— £) + &(— () ~ 
Va 2 + 6 2 (-0; a,b>0. Thus 

Z~N(o;T) ^ -£~AT(o ; r). 
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The following proposition and corollary show that A/"(0; T) is a uniquely de- 
fined sublinear distribution on (R n , Ci,Li P {R n ))- We will show that an A/"(0; T) 
distribution is characterized, or generated, by the following parabolic PDE 
defined on [0, oo) x R: 

Proposition 5.1.5 For a given tp G CiLi P (R d ) the function 
u(t,x) =R[(p(x + Vt£)], (t,x) G [0,oo) x R n 

satisfies 

u(t + s, x) = K[u(t, x + \/s£)] 

u(t, •) G CiLi P (R n ) for each t and u is locally Holder in t. Moreover u is the 
unique viscosity solution of the following PDE 

fin 

^-G(D 2 u)=0, u\ t=0 = <p, (5.2) 

here D 2 u is the Hessian matrix of u, i.e., D 2 u = {d 2 ixj u) d , =1 . 

Remark 5.1.6 In the case where £ has zero variance-uncertainty, i.e., there 
exists 70 G §(n) + 7 such that 

G{A)= l -{A, l0 ). 

The above PDE becomes a standard linear heat equation and thus, for G° = 
G{ 7 2} ; the corresponding G ^-distribution is just the d-dimensional classical 
normal distribution with zero mean and covariance 70 , i.e., £ ~ A/"(0,7o). 
In a typical case where 70 = Id G V we have 

%>(0] = tA^ / exp[-^^2]^(x)dx, V eC LLip (R d ). 

In the case where 70 G T from comparison theorem of PDE 

F^]>FaM Vy? G C LLzp (R d ). (5.3) 

More generally, for each subset fcT let 

G'(A) :=sup(A, 7 ) 
7Gr' 

the corresponding G-normal distribution A/"(0, T') is dominated by P G in the 
following sense: 

F Gr ,fa] -F Gr ,M <F Gr b-^], Vy?,^ G C LLrp (R d ). 
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Remark 5.1.7 In the previous chapters we have discussed l-dimensional 
case which corresponds to d = 1 and T = [a 2 ,^ 2 ] C R. In this case the 
nonlinear heat equation A5.2\) becomes 

du 

— - G{d 2 xx u) = 0, u\ t=0 = <pe C LLip {R), 

with G(a) = |(cf 2 a + — a 2 a~). In the multidimensional cases we also have the 
following typical nonlinear heat equation: 

i=\ 

where Gi(a) = |(cf 2 a + — £ 2 a _ ) and < a i < ~Ui are given constants. This 
corresponds to: 

T = {diag[j h ■ ■ ■ , 7J : 7, e [a- , of], i = 1, • • - , d}. 

x + Vt£ is J\f(x, ^-distributed. 

Definition 5.1.8 We denote 

P t G (<p)(x)=E[(p(x + VixZ)) = u{t,x), (£, x) € [0, 00) x R d . (5.4) 
Since for each (p 6 C^^IR ') we /iai>e £/ie Markov-Nisio chain rule: 

Remark 5.1.9 This chain rule was first established by Nisio LJ5\! and LJ^jj 
in terms of "envelope of Markovian semi-groups". See 155^ for an application 
to generates a sublinear expectation space. 

Lemma 5.1.10 Let £ ~ J\f(x,tT). Then for each a e M. d , the random vari- 
able (a, £) is A/"(0, [a 2 , a 2 ]) distributed, where 

a 2 = E[(a,0 2 ], a 2 = -E[-(a,0 2 ]. 

Proof. If £ is independent to £ and £ ~ £ then (a, £) is also independent to 

(a, C) and (a, () ~ (a, £). It follows from a£ + (3( ~ a/o 2 + /3 2 £ that 

a(a, + /3(a, C) = (a, «£ + [3Q ~ y/tf + 0*{&, £). 

Thus (a,O~AT(0, [a 2 ,a 2 ]). ■ 
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5.2 G— Brownian motions under G— expectations 

In the rest of this paper, we set Q = Cq(R + ), the space of all IR d -valued 
continuous paths (ou t ) teK +, with u>o = 0. Q is the classical canonical space 
and uu = (ut)t>o is the corresponding canonical process. It is well-known 
that in this canonical space there exists a Wiener measure (ft,!F,P) under 
which the canonical process B t (uj) = uot is a rf-dimensional Brownian motion. 
For each fixed T > we consider the following space of random variables: 

LUF T ):={X(w)=<p(wt 1 ,---,u>t m )- Vm>l,t ir --,t m e [0,7], 



V e C LLrp (R dxm )}- 

It is clear that {T~Ct}t>o constitute a family of sub-lattices such that 7i® C 
L° p (jF T ), for t < T < oo. L° p (jF t ) representing the past history of uj at the 
time t. It's completion will play the same role of Brownian filtration Tf as 
in classical stochastic analysis. We also denote 

00 
n=\ 

Remark 5.2.1 It is clear that C LLip (R d><m ) and then L q 1p (^t), L° p (F) are 
vector lattices. Moreover, since <p,ijj E Ci,Li P (R dxm ) imply ip-ip E Ci.Li P (R dxm ) 
thus X, Y E L° ip (F T ) imply X -Y E L? p (F T ); X, Y E L%{F) imply 
X-YEL%{T). 

We will consider the canonical space and set B t (uo) = oo t , t E [0, oo), for 
to E Q. 

Definition 5.2.2 The d- dimensional process (B t ) t > is called a G-Brownian 

motion under a sublinear expectation (Q,7i,K) defined on L® (T) if Bq = 
and 

(i) For each s,t>0 and ip E C LLip (R d ), B t ~ B t+S - B s ~ A/"(0, G); 
(ii) For each m = 1, 2, • • • , = t < t x < • • • < t m < oo the increment 
Bt m ~ B tm _ x is independent to B tl ,- ■ ■ , B tm _ ± . 

Definition 5.2.3 The related conditional expectation ofip(B tl —B to , ■ ■ ■ , B tm — 
Bt m _ 1 ) under 7i tk is defined by 

m[p{B tl - B t0 , ■■■ ,B tm - B tm )\H tk ] = <p m -k(Bt ir ■ • , B tk ), (5.5) 

where (p m ^ k {x l , ■ ■ ■ ,x k ) =E[ip(x 1 ,--- ,x k ,B tk+1 -B tk ,--- , B tm - B t J]. 
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It is proved in [55] that E[-] consistently defines a nonlinear expectation 
satisfying (a)-(d) on the vector lattice L^ p (J-t) as well as on L® (J 7 ). It 
follows that E[|X|], X G L? (Fr) (respectively L° p (F)) is a norm and thus 
L%{3~t) (respectively, L° p (jF)) can be extended, under this norm, to a Banach 
space. We denote this space by 

Ht = L^Ft), T G [0, oo), (respectively, H = L G (F)). 

For each < t < T < oo, we have L l G (F t ) C L^(F T ) C L^(F). In L G (T T ) 
(respectively, L g {!Ft))-, E[-] still satisfies (a)-(d) in Definition 13. 1.21 

Remark 5.2.4 i£ zs suggestive to denote L\ ' (Ft) fry 7Y° ; L G [Tt) by Tit and 
L G (J-) by 7i and thus consider the conditional expectation E[-|%] as a pro- 
jective mapping from 7i to Tit. The notation Lq^Ti) is due to the similarity 
with L 1 (n,Ft,P) in classical stochastic analysis. 

Definition 5.2.5 The expectation E[-] : L G (J-") i— > R introduced through the 
above procedure is called G— expectation, or G-Brownian expectation. The 
corresponding canonical process B is said to be a G-Brownian motion under 
E[-]. 

For a given p> 1 we also denote L P G (F) = {X 6 L G {F) : |X| P G Lq(F)}. 
L G {T) is also a Banach space under the norm ||X|| := (E[|X| p ]) 1/,p . We have 
(see Appendix) 

II V _i_ VII <r" II VII _i_ II VII 

||A +y\\ p < ||A|| p + yy || p 

and, for each X G L% Y G L^(Q) with ± + ± = 1, 

||XY||=E[|Xy|]<||X|| p ||X|| ? . 

With this we have ||X|| < ||X|| , if p < p'. 

We now consider the conditional expectation introduced in Definitior i5.2.3l 
(see ( 15. 5ft ). For each fixed t = tk < T the conditional expectation E[-|?Yt] : 
L° p (Fr) >— > I/° p (Ft) is a continuous mapping under ||-||. Indeed, we have 
E[E[X|W t ]] = E[X], X G £° P (F T ) and since if is subadditive 

E[X\H t ] - E[Y\H t ] < E[X - Y\H t ] < E[|X - F||ftt] 

We thus obtain 

E[E[X|^t] - E[Y|W t ]] < E[X - Y] 
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and 

t[X\Ht]-t[Y\Ht] < \\X-Y\\. 

It follows that E[-|W t ] can be also extended as a continuous mapping 

E[-\H t ] : L G {T T ) - L G {T t ) 

. If the above T is not fixed then we can obtain E[-|7i t ] : L G (!F) i— ► L^(jF t ). 

Proposition 5.2.6 We list the properties ofE[-\Tt t ], t G [0,T] t/iat hold in 
L° ip (F T ) and still hold for X,Y e L l G (F T ): 

(i) E[X\H t ] = X, for X G L G {F t ), t < T. 
(ii) IfX>Y then E[X\H t ] >E[Y\H t ]. 
(Hi) E[X\H t ] - E[Y\H t ] < E[X - Y\H t \. 

(iv) E[t[x\n t ]\n s ] = E[x\n tAs ], E[t[x\n t }} = e[x}. 

(v) E[X + v \H t ] = E[X\H t ] +r 1 ,r 1 e L G (F t ). 

(vi) E[r)X\H t ] = r] + E[X\H t ] + rj-E[-X\H t ] for bounded rj G L G (T t ). 

(vii) We have the following independence: 

E[X\H t ] = E[X], VX G L G (J^), VT > 0, 

where L G (J r ! r ) is the extension under ||-|| o/ L° p (jF^) which consists of ran- 
dom variables of the form (p(Bl , B\ 2 , • • ■ , B\ ), ip G CiL ip (M m ) , ti, • ' • ,t m G 
[0, T], m — 1, 2, • • • . Here we denote 

B l s = B t+s -B u s>0. 

(viii) The increments of B are identically distributed: 

E[<p(B$ i: Bl ■ ■ ■ , B$J] = E[<p(B tl: B t2: • • ■ , B t J). 

The meaning of the independence in (vii) is similar to the classical one: 

Definition 5.2.7 An R™ valued random variable Y G (L G (J-')) n is said to 
be independent to Tit for some given t if for each tp G Ci,Li P {^ n ) we have 

E[p(y)\7u] = t,[p(y)]. 
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It is seen that the above property (vii) also holds for the situation X G 
Lg-(jF') where Lq^F 1 ) is the completion of the sub-lattice Ut^o-^g^t) un der 

INI- 

From the above results we have 

Proposition 5.2.8 For each fixed t > (-B*) s >o is a G-Brownian motion 
in L G (JF') under the same G- expectation E[-]. 

Remark 5.2.9 We can also prove that the time scaling of B, i.e., B = 
(vAB t /\) t >o also constitutes a G-Brownian motion. 

The following property is very useful. 

Proposition 5.2.10 LetX,Y G L X G {T) be such that E[Y\H t ] = -E[-Y\H t ], 
for some t G [0, T\. Then we have 

E[X + Y\H t ] = fc[X\Ht] + W\H t ]. 

In particular, ifE[Y\H t ] = E[-Y\H t ] = then E[X + Y\H t ] = E[X\H t }. 

Proof. This follows from the two inequalities E[X + Y\H t ] < E[X\H t ] + 
E[Y\H t ] and 



E[X + Y\H t ] > E[X\H t ] - E[-Y\H t ] = E[X\H t } + E[Y\H 



t ■ 



Example 5.2.11 We have 

E[(AB t , B t )\ = a A t = 2G(A)t, VA G S(d). 
More general, for each s <t and rj = (r]^)fj =1 G L? G {T S \ S(d)) 

E[(r]B s t ,B s t )\n s ] = a v t = 2G{rj)t, s,t>0. (5.6) 

Definition 5.2.12 We will denote in the rest of this paper 

B*=(a,B t ), for each a = (oi, • • • ,a d ) T eR d . (5.7) 

Since (a.,B t ) is normal distributed (a.,B t ) ~ A/"(0, [cr aa r, c_ aa r]): 

E[tp(B?)} = E[^((a, Vim 

Thus, according to Definition \5.2.2\ for d-dimensional G-Brownian motion 
B a forms a 1 -dimensional G a Brownian motion for which the G a expectation 
coincides with EM. 
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Example 5.2.13 For each < s — t we have 

E[il>(B t -B s )\H a ]=fc[if>(B t -B s )]. 

If ip is a real convex function on K. and at least not growing too fast then 

E[X<p(B$ - Bf)\H t ] 
= X + E[<p(B* - B?)\H t ] + X-±[-<p(B$ - Bf)\H t ] 

X + f°° x 2 

( p(x)exp(——— )dx 



2vr(T - t)a aa r J-oo V 2(T - t)a aa r 

X' f°° x 2 

tp(x) exp( — - -)dx. 



v/27r(T-t)|a_ aa T|7-oo V ' 2{T-t)\<J_ 
In particular, for n — 1, 2, • 
E[\B?-BT\H S ]=E[\B?_ 



S\ J 

1 f°° X 2 

' x\ n exp(— — )dx. 



2n(t-s)a aaT J-oc 2(t-s)(T a 

But we have E[-|£f - Bf\ n \H s ] = E[-\B?_ s \ n ] which is when a_ aa r = 
and equals to 

1 f°° x 2 

l a;| n exp( -)dx, ifo-_ aa r<0. 



^27r(t-s)\a_ aaT \J^oo 2(t - s)\a_ aa T 

Exactly as in the classical cases we have E[Bf — Bf\7i s ] = and 

t[(Bf - Bf) 2 \H s ] = a aaT (t - s), E[(B? - Bf^Hs] = 3a 2 aaT (t - s) 2 , 
E[(B? - Bff\H s ] = 15 *l aT (t-s)\ E[(B? - Bf) 8 \H s ] = 105^ aT (t - S ) 4 , 

n\B? ~ Bf\\H s ] = vgg^SZ, E[\Bf - BffW ~ 2V ~ m ~ SVaaTf/2 



7T a/7T 



E\\Bf-B^\n,]=8y^- s)a - T]5/2 



7[ 



Example 5.2.14 For each n = 1, 2, • • • , < t < T and X € L^Ft), we 

have 

E[X{B* - B?)\H t ] = X + E[(B* - Bf)\H t ] + X~E[-(B^ - Bf)\H t ] = 0. 
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This, together with Proposition \5.2.10i yields 

E[Y + X(B^ - Bf)\Ht\ = E[Y\H t ], Y e L l G (F). 

We also have, 

E[X(B* - B?) 2 \H t ] = X + E[(B* - Bf) 2 \H t ] + X~E[-(B% - Bf) 2 \H t ] 
= [X + ( r SL3L T+X-a_ aa T](T-t). 

Remark 5.2.15 It is clear that we can define an expectation E°[-] on L\ ' (•T 7 ) 
in the same way as in Definition \5.2.2\ with the standard normal distribution 
JV(0, 1) in place o/JV(0, T). IfI d eT then it follows from (E3J that jV(0, 1) 
is dominated by Af(0,T): 

E°{v(B t )]-E°i4,(B t )]<Ei(v-?P)(B t )]. 

Then E°[-} can be continuously extended to Lq(F). E°[-] is a linear expecta- 
tion under which (B t )t>o behaves as a classical Brownian motion. We have 

- E[-X] < E°[X] < E[X], -E[-X\H t ] < E°[X\H t ] < E[X\H t ]. (5.8) 

More generally, ifT'cT since the corresponding P' = P Gl " is dominated by 
P G = P Gr , thus the corresponding expectation E' is well-defined in Lq(T) 
and E' is dominated by E: 

E'[X]-E'[Y]<E[X -Y], i,ye4(f). 

Such kind of extension through the above type of domination relations was 
discussed in details in [55]. With this domination we then can introduce a 
large kind of time consistent linear or nonlinear expectations and the corre- 
sponding conditional expectations, not necessarily to be positive homogeneous 
and/or subadditive, as continuous functionals in Lq^J 7 ). See Example \5.3.14 
for a further discussion. 

Example 5.2.16 Since 

E[2Bf(B? - Bf)\H s ] = E[-2Bf(B? - Bf)\H s ] = 

we have 

E[(B?y - m 2 \n s ] = t\(Bf - Bf + b:) 2 - m 2 m 

= E[(B? - Bf) 2 + 2(B? - Bf)Bf\H s ] 
= «r aa r(t-s) 
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and 

nmy - m 2 ) 2 \n s ] = ek(b? - b» + B*y - m^m 

= E[{(B? - Bf) 2 + 2(B? - Bf)Bf} 2 \H s ] 

= E[(B? - Bf)' + 4(5f - Bf) 3 B» + 4(£f - Bf) 2 (Bf) 2 \H s ] 

< t[(Bf - Bf)*} + 4E[\B* - Bf\ 3 ]\Bf\ + 4cr aa ,(t - s)(Bf) 2 



ZvLAt ~ s) 2 + 8^[a aaL r(t - s)f/ 2 \Bf\ + 4a aa T(t - s)(B 



f? 



5.3 Ito's integral of G— Brownian motion 

5.3.1 Bochner's integral 

Definition 5.3.1 For T G R+ a partition tit of [0, T] zs a finite ordered 
subset n — {ti, ■ ■ ■ , t^} such that = to < ti < • • ■ < t^ — T. 

}x(ir T ) = max{\t i+ i - U\ ■ i = 0,1,- ■ ■ ,N - 1}. 

We use Ti^ = {t^jt^,--- ,t%} to denote a sequence of partitions of [0, T] 
such that lim A r_ >00 /i(7Ty ) = 0. 

Let p > 1 be fixed. We consider the following type of simple processes: 
For a given partition {t , • • • , t^} = tit of [0, T] we set 

AT-l 

Vt(") = ^2^k(uj)l[t k ,t k+1 )(t), 

fc=0 

where £& G L v G [Tt k \ A; = 0, 1, 2, • • • ,N — 1 are given. The collection of these 
processes is denoted by Mq (0, T). 

Definition 5.3.2 For an 77 G Mg°(0,T) with Vt = Ek=o &Ml[t fcl t fc+ i) (*), 

£/ie related Bochner integral is 

N-l 



/ r)t(w)dt = ^2 Z.k(u)(t k +i - t k ). 
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Remark 5.3.3 We set for each 77 e Mg'°(0,T) 



1 /- T 1 JV_1 

'M : = 7f / E[7^]d< = - J2 E&Mfa+i - t k ). 

1 JO - 1 ,_ n 



It is easy to check that E^ : M G ' (0, T) 1 — > R forms a nonlinear expectation 
satisfying (a) -(d) of Definition \3.1~M. We then can introduce a natural norm 



MIJ. = Erlfol] = i y E[|,fc|]<ft. 



Under this norm M G (0,T) can 6e extended to M G (0,T) which is a Banach 
space. 



Definition 5.3.4 For each p > 1 we denote by M G (0,T) the completion of 
Mq (0, T) under the norm 

T / JV-1 \ 1 / p 

(~jf IIM P II^) 1/P =(^E^HI"K^-^) • 

We observe that, 

rn TV— 1 7 1 

n^l tk(»)dt\] < ^Yl II^MH (Wi-*») = ^/ E[|^|]dt. (5.9) 
We then have 



Proposition 5.3.5 T7ie linear mapping J r)t{uS)dt : M G ' (0, T) h-» L g (J-t) 
is continuous and thus can be continuously extended to M G (0,T) h-» Lq^t)- 
We still denote this extended mapping by J Q rj t (u)dt, r\ <G M^(0, T). 

Since Mq(0, T) C M^,(0,T) for p > 1, this definition makes sense for 
V eM G (0,T). 
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5.3.2 Ito's integral of G— Brownian motion 

We still use Bf : = (a, B t ) as in (157711 . 

Definition 5.3.6 For each rj G M G (0, T) mi/i t/ie /orm 

JV-l 

»7t(w) = ^^(^)I [tfc ,t fc+1 )W 

fc=0 

we define 

fT N-l 

Jfa) = / r}(s)dB» := J2 UB? k+1 - Bl k ). 
Jo k=0 

Lemma 5.3.7 We have, for each 77 G M G (0, T), 

E[[ V (s)dBf]=0, (5.10) 

Jo 

E[(/ r^W) 2 ] < <x aa r / % 2 ( S )]ds. (5.11) 

Jo Jo 

Consequently, the linear mapping I : M G (0,T) 1 — > L? g [Tt) is continuous 
and thus can be extended to I : M G (0,T) 1 — > L G (jF T ). 

Definition 5.3.8 We define, for a fixed r\ G M G (0,T) ; the stochastic calcu- 
lus 

[ V (s)dBf := I( V ). 
Jo 

It is clear that JOfl) and ( fOI]) sfcH hold for r\ G M G (0, T). 
Proof of Lemma 15.3.71 From Example 15.2. 141 for each k 

EMB? k+1 - B? k )\H tk } = 0. 
We have 

E[ / ^(a)dB*] = E[ / V (s)dB* + ^-1^ - B^.J] 
Jo Jo 



E[ / i;( a )dBJ + E[^_ 1 (Sf iv - B^JI^.J 
Jo 



E[ / i7(a)dflf]. 
Jo 
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We then can repeat this procedure to obtain (I5.10p . We now prove (15. lip 

E[(J V (s)dBf) ]=E[Q^ 1 V (s)dB: + ^ l (B: N -Bt N _ 1 ) ) j ] 

= E[ / V (s)dBl 



o 



+ E[2 Qf ^ T)(s)dB*\ 6v-i(£f„ - B? N J + e N -i(B? N - B^J'in^]] 
= E[(J N1 V (s)dBf\ + e N -i^T{t N - iiv-i)]. 

Thus E[( Jj w v(s)dBff] < E[ ( f^- 1 r / ( S )rf J Bf) 2 ] + E^.J^r^ - iiv-i). 
We then repeat this procedure to deduce 

E[(/ r ] (s)dB s ) 2 ]< ( T aaT J2mk) 2 ](t k+1 -t k )= [ E[(r/(t)) 2 ]dt 



fc=0 



We list some main properties of the Ito's integral of G-Brownian motion. 
We denote for < s < t < T 

t ,T 

VudB*:= / I {Stt] {u)ri u dB*. 

Js JO 

We have 

Proposition 5.3.9 Let rj, B e M%(Q,T) andO<s<r<t<T. Then in 
Lq(Ft) we have 

(i) JlnudBi = r; Vu dB- +jl Vu dB*. ^ 

(ii) J (an u + 9 u )dB* = a J rj u dB* + J 9 u dB*, if a is bounded and in Lq^J 7 ,,). 
(iii) t[X + £ TfrdBZlH.] = t[X\H s {,VX e L l G (F). 
(iv) E[(£ Vu dB*) 2 \Hs} < a aaT f*E[\r) u \ 2 \H s }du. 

5.3.3 Quadratic variation process of G— Brownian mo- 
tion 

We now consider the quadratic variation process of G-Brownian motion. 
It concentrically reflects the characteristic of the 'uncertainty' part of the 
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G-Brownian motion B. This makes a major difference from the classical 
Brownian motion. 

Let 71"^, N = 1, 2, • • • , be a sequence of partitions of [0, t}. We consider 



N-l 



to") 2 = £[(*&/ -(*#) 

fc=0 
N-l 



a \2 



■"fc+1 
fc=0 fc=0 



As £4(71"^) = max <fc<Ar_i(t^ fl — t^) — > 0, the first term of the right side 
tends to 2 J BfdBf. The second term must converge. We denote its limit 
by (B a ) v i.e., ' 



N-l r t 

( B ^t= I™ £(%, -B%f = {Blf-2 BfdBf. (5.12) 

By the above construction (B a ) t , t > 0, is an increasing process with (B a ) = 
0. We call it the quadratic variation process of the G-Brownian motion 
B a . Clearly (-B a ) is an increasing process. It is also clear that for each 
< s < t and each smooth real function ip such that ip((B a ) t _ s ) G L^(jF t _ s ) 

we have E{ip((B a ) t _ s )] = E[ip((B a ) t - (B a ) s )]. We also have 

(B a ) t = (B- a ) t = (-B a ) t . 

It is important to keep in mind that (B a ) t is not a deterministic process 
except in the case a aa r = —u_ aa T and thus B a becomes a classical Brownian 
motion. In fact we have 

Lemma 5.3.10 For each < s < t < 00 

E[(B a ) t - (B a ) s \H S ] = o aa r(t - s), (5.13) 

E[-((B a ) t - (B a ) s )\H s ] = a_ aa r(t - s). (5.14) 

Proof. From the definition of (B a ), Proposition I5.3.91 (iii) and Example 
15.2.161 

t[(B a ) t - (B a ) s \H S ] = t[(B a ) 2 - {B a s f - 2 f B a u dB a u \H s ] 

J s 

= t[(B a r-(B a ) 2 \H s ]=a aaT (t-s). 
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We then have ( 15.131) . ( 15.14p can be proved analogously by using the equality 
E[-((BfY-(Bff)\n s } = a^ aT (t-s). m 

An interesting new phenomenon of our G-Brownian motion is that its 
quadratic process (B) also has independent increments. In fact, we have 

Lemma 5.3.11 An increment of (B a ) is the quadratic variation of the cor- 
responding increment of B a , i.e., for each fixed s > 

(B a ) t+s -(B a ) s = ((BT) t , 
where B s t = B t+S - B s , t > and (B s )f = (a, B\). 

Proof. 

(B a ) t+S ~ (B a ) s = {B? +S f - 2 jf* BldBl - ((Bff - 2 J' B a dB* 

= (b? +s - B*y - 2 £ + \b: - B»)dB: 

= (b: +s - b:y - 2 [\b: +u - B*)d(B* +u - s») 

Jo 

= ((BT) t - 



Lemma 5.3.12 We have 

E[(B*) 2 } = E[((B a ) t+s - (B a )f\H s ] = ol &T t\ s,t>0. (5.15) 

Proof. We set <p(t) := E[(-B a ) t 2 ]. 

V (t) = t[{{Blf-2 f B*dB* u Y} 
Jo 

< 2E[(Bt) 4 } + 8E[( f BldBlf] 

Jo 

< 6a 2 aaT t 2 + 8a aa T / E[{B*) 2 ]du 

Jo 



lOo-Lrt 2 . 
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This also implies E[((£ a ) t - (B a ) s ) 2 ] = <p(t - s) < l0a 2 aaT (t - s) 2 . For each 

s e [o, t) 

V (t)=n((B a } s +(B*} t -(B*} s ) 2 } 

< t[((B*) s ) 2 ] + t[((B*) t - (B*) s ) 2 } + 2E[«P a ) t - <P a >J (B*) J 

= <p( 8 ) + <p(t -s) + 2E[E[«P a ) t - (B*) S )\H S ] (P a )J 
= <p(s) + <p(t-s) + 2cr 2 aaT s(t-s). 

We set 5n = t/N, t% = kt/N = k5^ for a positive integer N. From the 
above inequalities 



<p(t%) < ¥>(#-i) + <p{&n) + 2^ a -4-i^ 



A' 



< cp(!%_ 2 ) + 2<p(6 N ) + 2a 2 aaT {t N N _ l + t N N _ 2 )5 



N 



We then have 

<p(t) < N<p(8 N ) + 2a 2 aaT J2 tffor < Wt 2 a 2 aaT /N + 2a 2 aaT ^ $6 

fc=0 fc=0 

Let iV -> oo we have w(i) < 2ct 2 t f' sds = a 2 T t 2 . Thus E[(-B a )?l < a 2 T t 2 . 

' * ' £13. J U £13. L \ ' t J 33 

This, together with E^P 3 ) 2 ] > E°[(B a ) 2 } = v 2 aaT t 2 , implies ( KWf . In the 
last step, the classical normal distribution P°, or iV(0,7 7(f), 7o € I\ is 
chosen such that 



JV-l JV-l 

lN r 

JV- 



tr[7o7 aa ] = cr aaT = suptr[77 aa ]. 

7er 



Similarly we have 

E[((Z?*> t - (P a ) s ) 3 |K] = alr(t - s)\ (5.16) 

n((B a } t -(B a } s T\H s } = a aaT (t-sr. 

Proposition 5.3.13 Let < s < t, £ e L^tJF,) and X e L^J 7 ). Then 

E[X + £((B?) 2 - (P a ) 2 )] = ELY + {(ST ~ Bf) 2 ] 

= nX + a(B a ) t -(B a ) s )]. 
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Proof. From ( 15.121) and Proposition 15.2.101 we have 

E[X + £((5f ) 2 - (B*) 2 )} = E[X + £«£ a ) t - (B a ) s + 2 J* B^)} 

= nX + Z((B a ) t -(B a ) s )]. 

We also have 

E[X + d{{B?) 2 - (Bf) 2 )] = E[X + £{(Bf - B:) 2 + 2{Bl - B*)B*}] 

= E[X + Z(B?-Bff]. 



Example 5.3.14 We assume that in a financial market a stock price (S t )t>o 
is observed. Let B t = log(St), t > 0, be a 1- dimensional G-Brownian motion 
(d = 1) with T = [a*, a*}, with fixed a* G [0, |) and a* G [1, oo). Two traders 
a and b in a same bank are using their own statistic to price a contingent 
claim X = (B) T with maturity T . Suppose, for example, under the probability 
measure P a of a., B t (u)) t >o is a (classical) Brownian motion whereas under 
Pb of b, T}Bt(u)t>o is a Brownian motion, here P a (respectively, F^) is a 
classical probability measure with its linear expectation E a (respectively, Eb ) 
generated by the heat equation d t u = \d 2 x u (respectively, d t u = \d 2 x u). Since 
E a and Eb are both dominated by E in the sense of (3), they can be both well- 
defined as a linear bounded functional in Lq(J-). This framework cannot be 
provided by just using a classical probability space because it is known that 
(B) T = T, W a -a.s., and (B) T = ? ; Pb-a.s. Thus there is no probability 
measure on Q with respect to which P a and Pb are both absolutely continuous. 
Practically this sublinear expectation E provides a realistic tool of dynamic 
risk measure for a risk supervisor of the traders a and b: Given a risk position 
X G L l G (F T ) we always have E[-X\H t ] > E a [-X\H t ] VE h [-X\Ht] for the 
loss —X of this position. The meaning is that the supervisor uses a more 
sensitive risk measure. Clearly no linear expectation can play this role. The 
subset T represents the uncertainty of the volatility model of a risk regulator. 
The larger the subset T, the bigger the uncertainty, thus the stronger the 
corresponding E. 

It is worth considering to create a hierarchic and dynamic risk control system 
for a bank, or a banking system, in which the Chief Risk Officer (CRO) uses 
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E = E G for his risk measure and the Risk Officer of the ith division of the 
bank uses E J = E • for his one, where 

G{A) = - sup tr[^ T A], Gi{A) = - sup tr[^ T A\, r< c I\ % = 1, ■ ■ ■ , I. 

Thus E* is dominated by E for each i. For a large banking system we can 
even consider to create E JJ = E ij for its (i,j)th sub-division. The reason 
is: In general, a risk regulator's statistics and knowledge of a specific risk 
position X are less than a trader who is concretely involved in the business 
of the financial loss X . 

To define the integration of a process r] G Mq(0, T) with respect to d (B a ) 
we first define a mapping: 

T N-l 

QoAv) = / v(s)d(B*) s := I>« SS \ +1 -^ a >J : Mg°(^T) h* L\F t ). 

^° fc=0 

Lemma 5.3.15 For each rj e M^°(0,T) 

E[|Q ,T(r/)|]<a aa T / E[\Vs\}ds. (5.17) 

Jo 

Thus Qo t r '■ Mq (0, T) y—> L 1 (jF T ) is a continuous linear mapping. Con- 
sequently, Q ,t can be uniquely extended to M@(Q,T). We still denote this 
mapping by 

/ i7(s)d<B a >, = go,T(»7), veM^T). 

Jo 

We still have 

E[| f V (s)d(B a ) s \]<a aaT [ E[\Vs\]ds, Vr? e M^(0,T). (5.18) 

Jo Jo 

Proof. With the help of Lemma 15.3.101 (15.17ft can be checked as follows: 

iV-l N-l 

eii E^« sa W - ( B X^ ^ £%*i -^W - ( B X I***]] 

fc=0 A;=0 

JV-1 



£e[|&|]<w(*jm-i -t k 

k=0 

0"aaT / E[|jfr|]ds. 

Jo 
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■ 
We have the following isometry. 

Proposition 5.3.16 Let rj G M G (0,T). We have 

E[( T V (s)dBf) 2 ] =E[/% 2 ( S )d(5 a )J. (5.19) 

Jo Jo 

Proof. We first consider rj G M G ' (0, T) with the form 

JV-l 

Vt(u) = ^2tk{u)I[t k ,t k+1 ){t) 

k=0 

and thus j T V (s)dBf : = EEo'&W^ - Bf k ). By Proposition EZM we 
have 

ELY + 2UB? k+1 - B? k MBl i+i - BQ\ = ELY], for X G L&f), I ? k. 

Therefore 

r T /JV-l \ 2 JV-l 

E[( / ^W) 2 ] = E[ J] &(S» +i - B£) ] = MJ2 $(K +1 ~ B V 2 ]' 

J ° \k=0 J k=0 

This, together with Proposition 15.3.131 yields that 

rp fs[ — ^ rp 

mf v(s)dB:r]=nj2e k ((B% +i - (B*) tk )} = E[ [ ^ood^).]. 

./o fc=0 ' ./o 

Thus ( 15. 191) holds for 7? G M G ' (0, T). Then we can continuously extend this 
equality to the case 77 G M^(0, T) and obtain (15.191) . ■ 

5.3.4 Mutual variation processes for G— Brownian mo- 
tion 

Let a = (ai, • • • , ad) T and a = (oi, • ■ • , a^) T be two given vectors in M. d . We 
then have their quadratic variation processes (B a ) and (B a ). We then can 
define their mutual variation processes by 

{B^B-*) t := l -[{B* + B*) t -{B*-B-) t ] 



4 
1 



-[(B*+*) t -(B*-*)J. 
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Since (B a - a ) = (B^ a ) = (-B a - a ) we see that (B a ,B a ) t = (B a ,B a ) t . In 
particular we have {B a , B a ) = (B a ). Let nf, N — 1, 2, • • • , be a sequence of 
partitions of [0, £]. We observe that 

n-i 1 N-l 

E(^ +1 - B k^ B h ~ ty = i Ek^K - 5 ^) 2 - ( 5 C - 5 r) 2 ]- 

fc=0 fc=0 

Thus as h(tt^) ->0we have 

JV-1 

lim V (S3* - B% ) {B% - B% ) = (B a , £? a ) . 

N~*0 ^— ' fc+i fc fc+i fc X ,t 



k=0 

We also have 



[ B ^B a ) t = \\( B ^) t -( Ba -~*) t \ 



-[{Bf +a f -2 [ B a+a dB a+a - {Bf- a f + 2 [ Bf- a dB 
4 Jo Jo 



a-a,;ua-al 



= B a B a - / B a dB a - / BfdBf. 
Jo Jo 

Now for each rj e Mq(0,T) we can consistently define 

f\d(B^B' a ) s = \j\ s d{B a ^ a ) s -\j\ s d{B a -- a ) s . 
Lemma 5.3.17 Let r) N G M^'°(0, T), N = 1, 2, • • • , 6e o/7orm 

JV-1 



fc=0 



wift ^(tt^) — > and i] N ^ rj in Mq(0,T) as N — ► oo. T/ien we nave £/ie 
following convergence in Lq(J-'t): 

f T n-i 

- / T^d^i^. 
./o 
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For notational simplification we denote B l = B ei , the i-th coordinate of 
the G-Brownian motion B, under a given orthonormal basis (ei, • • • , e<f) of 
M. d . We denote by 

({B) t h = {B\Bi) t . 

(B) t , t > 0, is an §(n)-valued process. Since 



We have 

E[((B) t ,A)] = 2G(A)t, AeS(n). 

Now we set a function 

v(t, X) = E[ip(X + (B) t )], (t, X) G [0, oo) x S(rc). 

Proposition 5.3.18 v solves the following first order PDE: 

d t v = 2G(Dv), v\t=o = v?, 

where Dv = (9 x ..)f,-i- We also have 

v(t,X) =snpif(X + tA). 
Aer 

Proof. (Sketch) We have 

v(t + 5,X)= E[<p(X + (B) s + (B) t+8 - (B) 5 )] 
= t[v(t,X + (B) 5 )]. 

The rest part of the proof is as in one dimensional case. ■ 

Corollary 5.3.19 We have 

(B) t e tT := {t x 7 : 7 G T}. 

or, equivalently d t r((B) t ) = 0, where djj(x) = inf{|x — y\ : y G U}. 

Proof. Since 

E[d tr ({B) t )} = supd tr (tA) = 0. 
Aer 

It follows that d tr ((B) t ) = in L l G (F). ■ 
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5.3.5 Ito's formula for G— Brownian motion 

We have the corresponding Ito's formula of $(X t ) for a "G-Ito process" X. 
For simplification we only treat the case where the function $ is sufficiently 
regular. 

Lemma 5.3.20 Let $ G C 2 (M. n ) be bounded with bounded derivatives and 
{dl» x vQ}™ u=1 are uniformly Lipschitz. Let s G [0, T] be fixed and let X = 
{X 1 , • • • , X n ) T be an n-dimensional process on [s, T] of the form 

X? = X: + o"(t -s)+ r) vij ((B\ B') t - (B\ B') a ) + ^{B\ - B{\ 

where, for v = 1, • • • , n, i,j = 1, • • • , d, a u , rf^ and f3 vt i are bounded ele- 
ments of L|.(jF s ) and X s = (X%, ■ ■ ■ ,X™) T is a given M. n -vector in Lq^s). 
Then we have 

$(X t ) - $(X,) = [ d x ^(X u )(3^dBl + [ d x MX u )a u du (5.20) 



[d x ^(X u )rj^ + l -dl x ^{XuW i r]d (B\ Bi) u . 

Here we use the Einstein convention, i.e., the above repeated indices [i,v, i 
and j (but not k) imply the summation. 

Proof. For each positive integer N we set S — (t — s)/N and take the 
partition 

ir^ t] = {t%,t? r -- ,t%} = {s,s + S r -- ,s + N5 = t}. 

We have 

JV-l 



$(X t ) - $(X S ) = J>(* Cl ) - $(%)] 

fc=0 
JV-1 

= Y,[dx>*(X t s)(X$ r -xi 



k=0 

+ l[d 2 x ^$(X t »)(X? N -XfriXfo-XM+ri?]] (5.21) 

2 k fe+1 u k fc+1 k 



where 

Vk = [^$(^+0 fe (^ +i -x t ^ 
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with 4 G [0,1]. We have 

E[|rf|] = E[|[^$(X 4f + 4(% +i - Xq)) - ShrQiXq)] 

x (*£ +1 ~ x J)( x *r +1 - *£ )ll 

<ct[\X t , +i -X t M\ 3 }<C[5 3 + 5 3 /% 

where c is the Lipschitz constant of {9^ii X v ^}t,v=i- ^ n the ^ as ^ s ^ e P we use 
Example 15.2.131 and (15.161) . Thus ^ fc E[|^|] — > 0. The rest terms in the 
summation of the right side of (15.211) are £j + £ t with 

AT-l 

ef = £{a^(x tf )K(tf +1 - *f ) + i^((b*, s% - (b*, b%) 

fc=0 

+ /^(BJv - B^)} + ±8^$(X#)P d F*(I& -B'){B{ N -B{ N )} 

l fc+l l k 2 k k + 1 k k+1 k 



and 

1 ^ 

fc=0 

x [a v {t» +x - t») + V » lm ((B\ B m ) tN - (B\ B m ) tN )] 

l k+l L k 

+ K(Ci - ^) + ^ ((£*, ^V - (B\ Bi) tN W\B\„ - BL). 

l k+l l k k+1 "k 

We observe that, for each u G [t^, t^ +1 ) 

N-l 

E[\d x ^(X u ) -Y,d x ^(X t ,)I [t ^ t , +i) (u)\ 2 } 



k=0 

:'"li \- V 121 s~ n\x i x2l 

k 



E[\d x ^(X u ) - d x ^(X tN 



<c 2 E[\X u -X t M\ 2 ] <C[5 + 5 2 
J2dl x ^(X t „)I [t » it » • ,(•) -> ^,S(X) in M£(0,T). 



Thus ^f^S,,*^)!^ )(•) tends to ^m$(X) in Mg(0,T). Similarly, 



AT-l 



fc=0 
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Let iV — > oo. From Lemma f5. 3. 171 as well as the definitions of the integrations 
of dt, dB t and d (B) t the limit of £f in L G {Tt) is just the right hand side of 
(15.201) . By the next Remark we also have C,f — > in l? G {Ti). We then have 
proved (1QD1 . ■ 



Remark 5.3.21 In i/ie proo/ o/C/^ — ► in L G {Tt), we use the following es- 
timates: for^p N e M G '°(0,T) such that if;? = J2k=o f£l [t jv ,^ +1 )(t), andvr^ = 

{0 < t , • • • , t N = T} with lim^ jufar ) = and ^k=o %£D(**+i-*") < 
C, for all N = 1, 2, • • • , we have E[| ^J^, 1 £f (if +1 - tf ) 2 ] -► and, for any 
fixed a, a E.W 1 , 

JV-l JV-l 



k=0 

£%ffc(Ci-0 2 -o, 



A:=0 fe=0 

AT-l 



fc=0 



JV-l 

e[| E^« SS V +1 - ( si V)(Ci - tf)ll 

fc=0 
JV-l 

< £% fc "l(Ci - tf) -E[((i? a ) tf+i - (5 a ) tf )|?%]] 

fc=0 

JV-l 



E%£W(Ci-*n 2 -o, 

fc=0 



as well as 



JV-l JV-l 



e[| ^ef(Ci -^)(^ +1 -^ a r )i] < E^iefiKCi -tf)E[|B^ -s^|] 



fc=0 fc=0 

■JV-l 



^E%?IKCi-tf) 8/a -o 

7T z — ' 
fc=0 
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and 



k=0 
N-l 

k=0 
AT-l 

fc=0 
AT-l 



fc=0 



We now can claim our G-Ito's formula. Consider 

X? = XX+ foT,da+ f\^d{B\Bi) s + f (3»JdB 
Jo Jo Jo 



Proposition 5.3.22 Let a u , (3 U ^ and rf 1 ^ , v = 1, • • • ,n, i,j = 1, • • • , d be 

bounded processes of M^(0,T). Then for each t > and $ G L^ G {Tt) we 
have 

$(X t ) - $(X.) = [ d x »$(X u )^dBl + [ d x MXu)a u u du (5.22) 



+ / [d x ^{X u ) V ^ + \dl^X u )Pl i PJ]d{B\Bi) u 



Proof. We first consider the case where a, r] and (3 are step processes of the 
form 

JV-l 

%(w) = E^Hi [tfc , tfc+1 )W- 

fc=0 

From the above Lemma, it is clear that ( 15.221) holds true. Now let 

X? N = X»+ f a»fds+ [ V v s ij > N d(B\Bi) s + f ^' N dB{ 
Jo Jo Jo 
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where a , r] and (3 are uniformly bounded step processes that converge 
to a, 7] and j3 in Mq(0, T) as N — > oo. From Lemma [5.3.201 

$(Xf ) - $(X ) = f d x ^(X»){3^ N dBi + f d x MXu)< N du (5.23) 
Jo Jo 

Since 

E[|Xf^-Xf| 2 ] 

< C [ T {nW N - <) 2 ] + E[|itf" - <| 2 ] + E[((3* N - (3») 2 }}ds 
Jo 



We then can prove that, in Mq(0,T), 

c^X^R'* - <9^$(XK 

We then can pass limit in both sides of (15.231) to get (15.221) . ■ 
Example 5.3.23 

(B t , AS t ) = £ Ayfljfl* = 2 £[^ / BJdB^ + A l3 (B\ B*)J 

= 2YXAy f Bid& + (A,(B) t )}. 
Jo 



I. J 



Thus 

E[(A, (B) t )] = E[(B t , AB t )\ = 2G(A)t. 
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5.4 G— martingales, G— convexity and Jensen's 
inequality 

5.4.1 The notion of G— martingales 

We now give the notion of G-martingales: 

Definition 5.4.1 A process (M t )t>o is called aG— martingale (respectively, 
G-supermartingale, G-submartingale) if for each < s < t < oo, we 
have M t £ L^Tt) and 

t[M t \H s ] = M s , (resp., < M s , > M s ). 

It is clear that for a fixed X £ Lq(F) 'E[X\Ti t ] t >o is a G-martingale. In 
general how to characterize a G-martingale or a G-supermartingale is still 
a very interesting problem. But the following example gives an important 
characterization: 

Example 5.4.2 Let M el,^= (y? l )f=i e M^(0,T;R d ) andr/ = {rf j )j >j=1 G 
Mq(Q, T;S(d)) be given and let 

M t = M + [ VudBl+ [ rf£d{B\Bi) u - [ 2G( Vu )du, te [0,T]. 
Jo Jo Jo 

Then M is a G-martingale on [0, T] . To consider this it suffices to prove the 
case 7] e il4'°(0,T;§(cf)), i.e., 

N-l 

vt = ^vt k h k .t h+1 ){t)- 



k=0 

We have for s £ [tjv-i, ^jv], 

&[M t \H.] = M s + E[ V l_ i ((B\B^) t - (B\B*) s ) - 2G(r ]tN _ 1 )(t - s)\H s ] 

= M S + t[ V l_ t {Bl - B\){Bi - Bi)\H s ] - 2G( VtN _ 1 )(t - s) 

= M S . 

In the last step we apply the relation ( 15.61) . We then can repeat this procedure, 
step by step backwardly, to prove the result for any s £ [0, tjv-i]- 

Remark 5.4.3 It is worth mentioning that for a G-martingale, in general, 
—M is not a G-martingale. But in the above example when rj = then —M 
is still a G-martingale. This makes an essential difference of the dB part 
and the d (B) part of a G-martingale. 
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5.4.2 G— convexity and Jensen's inequality for 
G— exp ect at ion 

A very interesting question is whether the well-known Jensen's inequality still 
holds for G-expectation. In the framework of ^-expectation this problem was 
investigated in [4] in which a counterexample is given to indicate that, even 
for a linear function which is obviously convex, Jensen's inequality for g- 
expectation generally does not hold. Stimulated by this example jlD] proved 
that Jensen's inequality holds for any convex function under a (/-expectation 
if and only if the corresponding generating function g = g(t, z) is super- 
homogeneous in z. Here we will discuss this problem from a quite different 
point of view. We will define a new notion of convexity: 

Definition 5.4.4 A C 2 -function h : R i — ► R is called G-convex if the fol- 
lowing condition holds for each (y, z,i)6BxR d x S(d): 

G(h'(y)A + h"(y)zz T ) - h'(y)G(A) > 0, (5.24) 

where h! and h" denote the first and the second derivatives of h, respectively. 

It is clear that in the special situation where G(D 2 u) = \^u a G-convex 
function becomes is a convex function in the classical sense. 

Lemma 5.4.5 The following two conditions are equivalent: 

(i) The function h is G-convex. 

(ii) The following Jensen inequality holds: For each T > 0, 

E[h(<p(B T ))] > h(E[<p(B T )]), (5.25) 

for each C 2 -function ip such that h((p(Bx)) and (p(Bx) 6 Lq{Tt)- 

Proof, (i) ^> (ii): From the definition u(t,x) := P^Yp](x) = E,[ip(x + B t )] 
solves the nonlinear heat equation ( 15. 21) . Here we only consider the case 
where u is a C 1,2 -function. Otherwise we can use the language of viscosity 
solution. By simple calculation we have 

d t h{u{t,x)) = h'(u)d t u = ti(u(t,x))G{D 2 u(t,x)), 

or 

d t h{u{t,x)) - G(D 2 h(u(t,x))) - f(t,x) = 0, h(u(0,x)) = h{cp{x)), 
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where we denote 

/(*, x) = ti{u{t, x))G{D 2 u{t, x)) - G{D 2 h{u{t, x))). 

Since h is G-convex it follows that / < and thus h(u) is a G-subsolution. 
It follows from the maximum principle that h(P^(ip)(x)) < P^(h(ip))(x). In 
particular (15.251) holds. Thus we have (ii). 

(ii) =^(i): For a fixed (y, z,A) eRxR d x §(d) we set <p(x) := y + (x, z) + 
\ (Ax, x). From the definition of P t G we have d t (P^((f)(x))\t=o = G(D 2 (p)(x). 
With (ii) we have 

h(P t G (<p)(x)) < P?{h{ V )){x). 
Thus, for t > 0, 

- t [h(P t G ( V )(x)) - %>(*))] < \[P?{h{ V )){x) - %,(*))] 

We then let t tend to 0: 

h'( v (x))G(D 2 v (x)) < G(D 2 xx h( v (x))). 
Since D x {p(x) = z + Ax and D 2 x <p(x) = A we then set x = and obtain 

ra. ■ 

Proposition 5.4.6 TTie following two conditions are equivalent: 

(i) the function h is G -convex. 

(ii) The following Jensen inequality holds: 

E[h(X)\n t }>h(E[X\n t ]), £>0, (5.26) 

for each X <E L^F) such that h(X) G L^J 7 ). 

Proof. The part (ii) =^(i) is already provided by the above lemma. We can 
also apply this lemma to prove (15.261) for the case X e L^J 7 ) of the form 
X = <p(B tl , • ■ ■ , B tm — B tml ) by using the procedure of the definition of E[-] 
and E[-|7Yt] given in Definitions 15.2.21 and 15.231 respectively. We then can 
extend this Jensen's inequality, under the norm ||-|| = E[| • |] to the general 
situation. ■ 
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Remark 5.4.7 The above notion of G -convexity can be also applied to the 
case where the nonlinear heat equation A5.2fy has a more general form: 

fin i 

— -G(u,Vu,D 2 u) = 0, u(0,x)=ip{x) (5.27) 

(see Examples 4-3, 4-4 an d 4-5 in JEE] )- In this case a C 2 -function h : 
R i — > R is said to be G -convex if the following condition holds for each 
(y,z,A) gRx R d x S(d).- 

G(y, h'(y)z, h\y)A + h"(y)zz T ) - ti(y)G(y, z, A) > 0. 

We don't need the subadditivity and/or positive homogeneity of G(y,z,A). 
A particularly interesting situation is the case of g-expectation for a given 
generating function g = g(y,z), (y,z) G R x M. d , in this case we have the 
following g -convexity: 

\h"{y)\z\ 2 + g(h(y), h'(y)z) - h'(y)g(y, z) > 0. (5.28) 

This situation is systematically studied in Jia and Peng [31]. 

Example 5.4.8 Let h be a G -convex function and X e L^(jF) such that 
h(X) e L l G (F). Then Y t = hCk[X\H t \), t > 0, is a G-submarUngale: For 
each s <t, 

t[Y t \H s ] = E[hCE[X\F t ])\F 8 ] > h(t[X\F s ]) = Y s . 



5.5 Stochastic differential equations 

We consider the following SDE driven by G-Brownian motion. 

X t = X + [ b(X s )ds+ [ h il (X s )d(B\B j ) s + [ a j (X s )dBl te[0,T], 
Jo Jo Jo 

(5.29) 

where the initial condition Xq e R™ is given and 

b,hij,aj : R n ^R n 

are given Lipschitz functions, i.e., \(p(x) — t fi(x')\ < K\x — x'\, for each x, 
x' G R n , if = b, rjij and o~j, respectively. Here the horizon [0,T] can be 
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arbitrarily large. The solution is a process X e M^.(0,T;R n ) satisfying the 
above SDE. We first introduce the following mapping on a fixed interval 
[0,T]: 

A.(T) := Y E il4(0, T; R n ) i — > M&(0, T; R") 

by setting A t = X t , t G [0, T], with 

A t (r)=X + X + / b(Y s )ds + [ h l3 (Y s )d(B\Bi) s + [ a J (Y s )dBl 
Jo Jo Jo 

We immediately have 

Lemma 5.5.1 For each Y,Y' e Mg(0,T;IR n ) we have the following esti- 
mate: 

E[|A t (Y) - A t (r)| 2 ] < c / E[|y s - y;i 2 ]rf S) t e [o, r], 

where the constant C depends only on K , V and the dimension n. 

Proof. This is a direct consequence of the inequalities (15.9p . (15.111) and 

(E3ED. ■ 

We now prove that SDE (15.291) has a unique solution. By multiplying 
e -2Ct on ^q^^ sides f the above inequality and then integrating them on 
[0,T], it follows that 

/ E[\A t (Y) - A t (Y')\ 2 ] e - 2Ct dt < C f e- 2Ct I t[\Y s -Y' s \ 2 }dsdt 
Jo Jo Jo 

= C f f e- 2Ct dtt[\Y s -Y^]ds 



(2C)- 1 C [ (e- 2Cs -e- 2CT )E[\Y s -Y;\ 2 ]ds. 
Jo 



We then have 

/ E[\A t (Y)-A t (Y')\ 2 }e- 2Ct dt<]- [ E[\Y t - Yl\ 2 ]e~ 2Ct dt. 
Jo * Jo 

We observe that the following two norms are equivalent in Mq(0, T; MJ 1 ) 



[ E[\Y t \ 2 ]dt~ [ E[\Y t \ 2 }e- 2Ct dt. 
Jo Jo 



From this estimate we can obtain that A(Y) is a contraction mapping. Con- 
sequently, we have 
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Theorem 5.5.2 There exists a unique solution X e Mg.(0,T;M") of the 
stochastic differential equation Ii5.29\) . 
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Chapter 6 

Other Topics and Applications 

6.1 Nonlinear Feynman-Kac formula 

Consider SDE: 

dX f s > x = b{Xl> x )ds + h(Xl> x )d (B) s + a(X l f)dB s , s e [t, T], 
X l f = x. 



Y? = E[*(X*f) + I g(X t r >*,Y r t ' x )dr+ I f(X? X*)d (B) r \H,\. 

J s J s 

It is clear that u(t,x) := Y t ' x G T{q, thus it is a deterministic function of 
(t,x). u(t,x) solves 

d t u + sup{(6(:r) + h(x)j, Vw) + \ (o~(x)^cr T (x), D 2 u) 
+g(x,u) +/(x,w)7} = 0, 

u\t= T = $. 

Example 6.1.1 Let B = (B 1 ^ 2 ) be a 2-dimensional G-Browian motion 
with 

G(A)=G 1 (a 11 ) + G 2 (a 22 ) 

Gi(a) = -ia 2 a + - a 2 a~) 
In this case by Ito's formula 

99 
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dX s = fiX a ds + vX s d (B 1 ) s + aX s dB 2 s , X t = x. 



u(t,x) = E[<p(X$*)] = E[p{X^)\K 



We have 



u(t,x) = E[u(t + 5,X t t '* s )]. 
From which it is easy to prove that 

d t u + sup (n + wy)xd x u H sup [d xx u] = 0. 

6.2 Markov- Nisio process 

Definition 6.2.1 A n- dimensional process (X t ) t > on (Q,TC,(TCt)t>OjE) is 
called a Markov-Nisio process if for each < s < t and each tp e C^iW 1 ) 
there exists a ip E Cb{M. n ) such that 

E[<p(X t )\H.]=i/>{X 8 ). 

Remark 6.2.2 When E is a linear expectation then this notion describes a 
classical Markovian process related to Markovian group. Nisio 's semigroup 
(see Nision [j^i , 14®! ) extended Markovian group to sublinear cases to de- 
scribe the value function of an optimal control system. 

6.2.1 Martingale problem 

For a given function F(x,p,A) : M. d x M. d x S(d) i — ► K. we call a sublinear 
expectation space solves the martingale problem related to F if there exists 
a time consistent nonlinear expectation space (Q,Tt,(Ti t )t>o,E) such that, 
for each cp e C£°(R n ) we have 

M t := cp(X t ) - ip{Xo) - I F(X S , D<p(X a ), D 2 ^(X s ))ds, t > 

Jo 

is a E-martingale. (see P. [SI], [55]). 
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6.3 Pathwise properties of G-Brownian mo- 
tion 

It is proved that 

E[X] = supE Q [X] 
QeP 

where V is a family of probabilities on (fl, B(Q)), Q = C o (0, oo; M. d ). We set 

c{A) = sup Eq[1 a ] 
Q&P 

c(-) is a Choquet capacity. We have 

Theorem 6.3.1 (Denis-Hu-Peng) There exists a continuous version of G- 
Brownian motion: We can find a pathwise process (B t (u))t>o on some Q C Q, 
with c(Q c ) = 0, such that, for each u G Q, B t {u) G Co(0, oo) and E,[\B t — 
B t \] = 0, Vte[0,oo). 
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Chapter 7 
Appendix 



We will use the following known result in viscosity solution theory (see The- 
orem 8.3 of Crandall Ishii and Lions [TT]). 

Theorem 7.0.2 Let m <EUSC((0,T) x Q t ) for i = 1, • • • ,k where Qi is a 
locally compact subset of M. Ni . Let ip be defined on an open neighborhood 
of (0, T) x Q\ X • • • X Qk and such that (t, xi, • • ■ , x k ) is once continuously 
differentiable in t and twice continuously differentiate in (xi, • • • , Xk) GQiX 
• • • x Q k . Suppose that t G (0, T), Xi G Qi for i = 1, • • • , k and 



w(t, xi, ■ ■ ■ , x k ) := ui(t, Xx)-\ h u k (t, x k ) - ip(t, xi, 

<w(i,x 1 ,--- ,x k ) 



Xk) 



for t G (0,T) and X{ G Qi. Assume, moreover that there is an r > such 
that for every M > there is a C such that for i — 1, • • • , k 



bi < C whenever (bi,qi,Xi) G V 2 ' + Ui(t,Xi), 

\xi — Xi\ + \t — t\ < r and \ui(t,Xi)\ + \qi\ + \\Xi\\ < M. 

Then for each e > 0, there are Xi G S(iVj) such that 

(i) (bi, D x .ip(t, xi, • • • , x k ),Xi) e V ' Ui(t,Xi), i = 1, • • • , k; 

(H) 

" X x ■ ■ ■ 



(7.1) 



-(- + PII)< 







x k 



<A + eA 2 



(Hi) bi-\ h b k = <f t (i, xi, • ■ ■ , Xk) 

where A = D 2 <p(x) G S kN . 
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Observe that the above conditions ( 17.11) will be guaranteed by having m be 
subsolutions of parabolic equations given in the following theorem. 

Theorem 7.0.3 (Domination Theorem) Let m-order polynomial growth 
functions Ui €.USC([Q,T] x M. N ) be subsolutions of 

d t u - Gi{D 2 u) = 0, z = l,...,Jfc, (7.2) 

on (0, T) x M. N . We assume that {Gi}f =1 satisfies the following domination 
condition: 

k k 

Y^ G ( X i) ^ ' f° T dl X i G § ( iV )' SU0h that Xl Xi - ' ( 7 ' 3 ) 

i=l j=l 

and G\ is monotonic, i.e., G\{X) > G\{Y) if X > Y. Then the following 
domination holds: If the initial condition satisfies wi(0, x) + ■ - • Wfe(0, x) < 
for each x G ~R N then we have 



«i(t, x) + • • • u k (t, x) < 0, V(i, x) G (0, T) x 



dA j 



Proof of Theorem 17.0.31 We first observe that for 5 > 0, since D 2 \x\ 2m > 
0, the function defined by U\ := ui — 5/(T — t) — 5\x\ 2m is also a subsolution 
of ( 17.21) with £ = 1, and satisfies with a strictly inequality: 



d t «i - Gi(£rui) = ft«i - Gi(D\ + &£>> 



(T-t) 



2 



< d t u x - G 1 (D 2 u 1 ) - — J— - < 



(T-t) 2 - (T-£) 2 ' 

Since Mi + U2 + ■••«&< follows from £ti + «2 + • • • Uk < in the limit 6 | 0, 
it suffices to prove the theorem under the additional assumptions 

d t m - G^D 2 ^) < -c, c := 5/T 2 and 

\im t ^T,\x\^ooUi(t, x) = — oo uniformly in [0,T) x R . 

To prove our result, we assume to the contrary that wi(s, z) H h«fe(s, z) = 

5 > for some (s, z) G (0, T) x R^ and 5 > 0. We will apply Theorem [7702] 

for x = (xi, • ■ ■ , Xk), Xi G M N and 

k k-1 

w(t,x) := ^Mi(t,Xj), (p a (x) := ~(^2\x i+1 -Xi\ 2 + |x fc -Xi| 2 ). 

i=l i=l 
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Since for each large a > the maximum of w — <p a achieved at some (t a , x a ) 
uniformly inside a compact subset of [0, T) x M fcxiV . Set 

k 
i=l 

It is clear that M a > 5. We can also check that (see [17] Lemma 3.1) 

(i) lim a ^ OQ (p a (t a ,x a ) = 0. 

(ii) lim^oo M a = linia^o «i(t a , xj) + ■ ■ ■ + u k (t a , x%)) (7.4) 

= sup (4vr [wi(t, x)-\ h u k (i, x)]. 

where (t,x) is a limit point of (t a ,Xi). 
If t a = 0, we have 



<5 < M a < J2ui(0,xf)-<f a (t a ,x a ). 



i=i 



But by J2i=i u i{^i x t) ~~ * Ylii=i u i^^i) — 0- So t a must be strictly positive 
for large a. It follows from Theorem 17.0.21 that, for each e > there exists 
X { e S(N) such that 

(b?,D x Mt a ,x a ),X t ) e J3?Ui(t a ,x?) for i = I,--- ,fc. 

and such that ^ i=1 6f = 0, 

/ X 1 ... \ 



um < 







. X fe _! 

. o x k I 



< A + eA 2 



where A = D 2 ip a (x a ) G S 3N is explicitely given by 

'In ■ ■ ■ —In —In \ 
A = aJ 3N + oilw, J3N = 



-In 
\ -In 



In —In 
-In In ) 



(7.5) 



106 CHAPTER 7. APPENDIX 

The second inequality of ( 17.51) implies ^2 i=1 Xi < 0. Setting 
pi = D Xl (p a (x a ) = a(2x" -x$- x%), 

p k = D Xk ip a (x a ) = a(2xl - xl_ x - x"), 
We have J2i=iPi = ° 5 Ya=i h = and (b^p^Xi) G Jg + Mi(t a , xf). Thus 

b l -G 1 (X 1 )<-c, 
bi-Gi(Xi)<0, z = 2,3,... ,k. 

This, together with the the domination condition ( 17. 3D of Gj, implies 

fc fc 

_ c = _ ^ _ c > _ ^ G .( Xj ) > 0. 

2=1 1=1 

This induces a contradiction. The proof is complete. ■ 

Corollary 7.0.4 (Comparison Theorem) Let G±,G 2 '■ S(N) i— ► R 6e given 
functions bounded by continuous functions on S(iV) and let G\ be monotone 
and 

Gi(X) > G 2 (X), VX G S(iV). 

Let «2 ELSC((0,T) x IR^) 6e a viscosity supersolution of d t u — G\(D 2 u) = 
and u 2 £USC((0,T) x M. N ) be a viscosity subsolution of d t u — G 2 (D 2 u) = 0, 
such thatui(0, x) > w 2 (0, x), for allx G M. N . Then we have u\{t,x) < U2{t,x) 
for all (t,x) G [0,oo) x R N . 

Proof. It suffices to observe that —U\ GUSC((0,T) x M. n ) is viscosity sub- 
solutions of d t u — G*(D 2 u) = 0, with G*(X) := —Gi(—X). We also have, 
for each X x + X 2 < 0, 

G,(X 1 ) + G 2 (X 2 ) = G 2 (X 2 ) - G 1 (-X 1 ) 
< G 1 (X 2 ) - G 1 (-X l ) 
<Gi(Xx + X 2 ) <0. 

We thus can apply the above domination theorem to get —u\ + u 2 < 0. ■ 
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Corollary 7.0.5 (Domination Theorem) Let Gi : §(JV) i— ► R ; i — 0, 1, be 
two given mappings bounded by some continuous functions on S(iV) and /et 
■Uj ELSC((0,T) xR^) be viscosity supers olutions of d t u — Gi(D 2 u) = respec- 
tively for % — 0, 1 and /e£ W2 ^e a viscosity subsolution of ' dtu — G\(D 2 u) = 0. 
VFe assume that G zs monotone and that G dominates G\ in the following 
sense: 

G l (X)-G l (Y)<G (X-Y), VX,Ye$(N). 

Then the following domination holds: If 

u 2 (0,x)-u 1 (0,x) <u o (0,x), Vi £ R N , 

then u 2 (t, x) — u±(t, x) < u (t, x) for all (t, x) G [0, oo) x ~R N . 

Proof. It suffices to observe that —U{ GUSC((0,T) x M n ) are viscosity 
subsolutions of d t u — G\(D 2 u) = 0, with G\(X) := —Gi(—X). We also have, 
for each X 1 + X 2 + X 3 < 0, 

G%X 1 ) + Gl(X 2 ) + G\X 3 ) = G 1 (X 3 ) - G 1 (-X 2 ) - G (-X 1 ) 

<G (X 3 + X 2 )-Go(-X 1 )<0. 

We thus can apply the above domination theorem to get u 2 (t, x) — Ui(t, x) < 
u (t,x). m 
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